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ABSTRACT 

We present the first central velocity dispersion (ao) measured from the 0.85/xm Calcium II Triplet 
(CaT) for 8 advanced (i.e. single nuclei) local (z < 0.15) Ultraluminous Infrared Galaxies (ULIRGs). 
First, these measurements are used to test the prediction that the "cr-Discrepancy," in which the 
CaT CTo is systematically larger than the CTq obtained from the 1.6 or 2.3/ini stellar CO band-heads, 
extends to ULIRG luminosities. Next, we combine the CaT data with rest-frame /-band photometry 
obtained from archival Hubble Space Telescope data and the Sloan Digital Sky Survey (SDSS) to 
derive dynamical properties for the 8 ULIRGs. These are then compared to the dynamical properties 
of 9,255 elliptical galaxies from the SDSS within the same redshift volume and of a relatively nearby 
(z < 0.4) sample of 53 QSO host galaxies. A comparison is also made between the /-band and //-band 
dynamical properties of the ULIRGs. We find four key results: 1) the tr-Discrepancy extends to ULIRG 
luminosities; 2) at /-band ULIRGs lie on the Fundamental Plane (FP) in a region consistent with 
the most massive elliptical galaxies and not low-intermediate mass ellipticals as previously reported 
in the near-infrared; 3) the /-band M/L of ULIRGs are consistent with an old stellar population, 
while at //-band ULIRGs appear significantly younger and less massive; and 4) we derive an /-band 
Kormendy Relation from the SDSS ellipticals and demonstrate that ULIRGs and QSO host galaxies 
are dynamically similar. 

Keywords: galaxies: evolution — galaxies: formation — galaxies: interactions — galaxies: peculiar — 
galaxies: kinematics and dynamics — quasars:general 



1. INTRODUCTION 

\ Questions about the formation and evolution of galax- 
ies are as challenging today as when the 3rd Earl of 
Rosse first sketche d his observations of "external neb- 
ulae" (jRossd I185O0 . Referring to M51, he remarked 
that their complexity and striking beauty could hardly 
be the result of static processes. Objects like M51 and 
the Antennae (NGC 4038/39) have been the focus of as- 
tronomical investigations since their appearanc e in Her- 

■ schel's Catalogues of Nebulae and Clusters ()Herschell 
' I1786D . As photographic plates replaced pencil and pa- 
per, surveys continued to catalog peculiar "external neb- 
ulae" with everincreasin^^ 

; Ce.g. IPeasdllQlTl [19201: Wer^\wBL IRedman fc ShSml 
\ |1938() . Observational work, including morphological clas- 
sification and measuremen t of dynamical properties (e.g. 

■ |Hubble''1930; ■Zwickv"1956D , the first A^-body simulations 
(e.g. iHplmberg 1941) and subseq uent numerical simula- 
tions (iAlTadinlll965l : IWrightil972t ) explored the possibility 

^ National Research Council Postdoctoral Fellow at Naval Re- 
search Laboratory, Code 7211, 4555 Overlook Ave SW, Wash- 
ington D.C. 20375, USA 

Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218, USA 

^ George Mason University, Department of Physics &: Astron- 
omy, MS 3F3, 4400 University Drive, Fairfax, VA 22030, USA 

^ Leibniz-Institut fiir Astrophysik Potsdam (AIP), An der 
Sternwarte 16, 14482, Potsdam, Germany, brothberg@aip.de, 
dr.barry.rothberg@gmail.com, current address 

^ Naval Research Laboratory, Code 7211, 4555 Overlook Ave 
SW, Washington D.C. 20375, USA 

^ European Southern Observatory, Alonso de Cordova 3107, 
Casilla 19001, Vitacura, Santiago, Chile 

^ Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, 
HI 96822, USA 



that peculiar galaxies represented the transformation of 
galaxies from one form into another by means of inter- 
action. 

These earlier works all led directly to the Toom re Hy- 
pothesis (jToomre &: ToomrelfT97l iToomrd [19771 ) which 
posits that when gas-rich spirals collide and merge to- 
gether they form a new, more massive elliptical galaxy 
and that this process is responsible for the formation 
of all or most ellipticals in the Universe. The gravita- 
tional interaction between the two spirals rearranges the 
stellar orbit s from circular to random via violent relax - 
ation (e.g., ILvnden-Belll [19671: IHiorth fc Madsea[T99l . 
The process of gaseous dissipation funnels gas into the 
common gravitational center of the coalescing system, 
which triggers intense star-formation deep within molec- 
ular cloud s and adds substantial mass to the final rem- 
nant (e.g..lBarnes fc Hernauistill99ll : lMihos fc Hernguisti 
ll994alll996D . The most intense mergers are Ultralumi- 
nous Infrared Galaxies (ULIR Gs), systems with Ltr (8- 
1000 Atin) > 10^^ (e.g. see iSanders fc Mirabell[l996l : 
ISandersI [19991 : 1.ToseDhl[l999l for a review). The hypoth- 
esis that there is a natural evolution from ULIRG to 
QSO is based on the idea that gaseous dissipation fuels 
more than a nuclear starburst. The accretion of both 
gas clouds and stellar remnants fu els the formation of 
an active galactic nucleus or AGN ([Sanders et al.|[r988l 
hereafter S88). This was further supported by the simi- 
larity between the observed bolomctric (Lboi) luminosi- 
ties and space- densities of ULIRGs and QSOs out to 
at le ast z~0.4 ([Soifer et al.lll986l: iCanahzo fc StocktonI 
120011) . Only after the obscuring medium in the ULIRG 
is cleared by (presumably) radiation pressure and su- 
pernovae explosions does the QSO become visible (e.g. 
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iHopkins et al|[2006l ). 

Photometric observations of spiral-spiral mergers, in- 
cluding ULIRGs and their lower luminosity (10^^ Lq 
< LiR < 10^^ Lq) counterparts, Luminous Infrared 
Galaxies (LIRGs), demonstrated strong evidence sup- 
porting the Toomre Hypothesis. This includes confir- 
mation from optical to near-IR wavelengths that the 
global surface brightness (SB) profiles of advanced merg- 
ers foll ow the same de Vau couleurs r^/'* stellar light 
profile (|de Vaucoul eurs' lQSSl that characterizes ellipti- 
cal galaxies (e.g., iSc hwcizer .198^ ^ Wright_ct al. 1990; 
Lutzj 119911: iSchweizeil IIQM Iffibbard fc van GOTkMnl 



1991 iVeilleux et al.l l200T"rRothberg fc JosepH |2004t 



Veilleux et all l2006t iRothberg fc Fischer! I2010bl here- 
after Paper I), i n line w ith predictions from numerical 
simulations (e.g., [Barnes 1988, 1992). Moreover, numer- 
ical simulations show that gaseous dissipation during the 
merger will form a rotating gas disk which undergoes a 
strong st arburst and transforms in t o a rotating stellar 
disk (e .g. E arne s fc Hernqui stlll991l: iMihos & HernquistI 
1994bt iBarnes fc Hcrnauis tril996t [jvlihos & Hcrnquist 

20081) . The star- 



19961: lBarnesrT2002t iHopkins et al.l 



burst generates a luminosity spike at small radii, r < 
1-2 kilopars ecs (kpc), in the surface b r ightness profile s 
of mergers (jMihos fc HernquistI [T994at iSpringell [2000l ). 
This excess light was first directly detected in the K- 
band SB profiles of mergers, including (U)LIRG s, and 
found to have Lk ~ io9-5 -io.5 f Rpthberg fc Joseph! 
!2004 hereafter RJ04). 'Hopkins etall (|2008l) mod- 
eled the same sample and demonstrated that the ex- 
cess light from younger stars alone could account for 
30% of the total stellar mass. (jHaan eFall [20T1I) 
detected similar excess light at iJ-band using HST 
for a larger sample of (U)LIRGs spanning a wider 
range of merger stages. Their Figure 14 appears to 
show an evolution of the excess light as a function 
of the merger stage including what could be peaks 
at first passage and final coalescence. These proper- 
ties, taken together with the observed vast quantities 
j^qQ-io Ttf^) of cold mol e cular g as (e.g . Solomon et afl 



T99I 'Brvant fc ScoviUc' TM^; "Solomon ct al. 1997 



Scoville et al,.1997.: Downes fc Solomon 1998: lono et aP 
2005! : !Greye et al.! 2009!) and vigorous star-formation 



rates (e.g.. lPrestwich et al.lll99^ lAnantharamaiah et al.! 
[2OOOI SFR), make ULIRGs prime candidates for the pro- 
genitors of giant ellipticals (gEs) and QSO hosts. 

However, a significant challenge arose for the Toomre 
Hypothesis and the S88 scenario when dynamical 
masses (Moyn) obtained from central velocity disper- 
sions ((To) using the 1.6 or 2.3 /im CO band-heads 
(hereafter denoted as (To.co) Siiid imaging at H (1.6 
fiui) or if-band (2.2 fj,m) implied that (U)LIRGs 
were the progenit ors of low to intermediate mass el- 
lipticals (e.g . iShier etHI 119961: iShier fc Fischeil 119981: 
James et alT 11999!: iGenzel et al.l 120011: iTacconi et al.! 



2002! : iDasvra et al.!l2006D . The near-infrared was used 



because it is less af fected by the p r esenc e of dust than op- 
tical wavelengths. iPu nlop et al.! (|2003l ) (hereafter D03) 
compared near-IR photometry, by assuming a fiducial 
(R—K) transformation, with optical imaging of nearby (z 
< 0.4) radio loud and radio quiet QSOs (RLQ and RQQ, 



respectively j3- D03 concluded that the small half-light 
radii of ULIRGs precluded them as candidates f or the 
progenitors of QSO host galaxies. IVeilleux et al.l (|2006D 
made a similar comparison using Hubble Space Telescope 
(HST) imaging obtained with the Near Infrared Camera 
and Muhi-Object Spectrometer (NICMOS) at F160W 
(~ fl-band) of ULIRGs and Palomar Green (PC) QSOs, 
along with QSOs from D03 (transformed from i?-band 
to F160W), and relatively nearby ellipticals. They con- 
cluded that ULIRGs were the progenitors of 1-2 L* ellip- 
ticals or SOs. 

Interestingly, when the Calcium II Triplet stellar ab- 
sorption lines (A ~ 0.85 ^m) were used to measure 
(Jo (hereafter (To,CaT) a very different picture emerged. 
A comparison between (To,caT and ao,co in the same 
set of LIRGs systernatical ly showed (Jo.CaT > <^o,co 
(IRothberg fc Joseph! [2006al hereafter RJ06a). More- 
over, the CTo.CaT values were consistent with ellipticals 
ove r a large mass range, including gEs. Rec ent work 
by (|Rothberg|l20^ IRothberg fc Fischei!l2010bl hereafter 
Paper I) has effectively explained this a -Discrepancy. 
Paper I compared (To.CaT with cto.co and /-band pho- 
tometry with if-band photometry in advanced LIRG and 
non-LIRG mergers as well as elliptical galaxies. No a- 
Discrepancy was found f or elliptical galaxies, a result 
subsequentlv confirmed bv lVanderbeke et al.! (|201l! ). Pa- 
per I also showed that in advanced mergers (Log Lir < 
11.99), the (T-Discrepancy strongly correlated with Log 
LiR and dust mass. Although Paper I did not include 
any ULIRGs in the sample, the results were extrapolated 
to brighter Ljr. Paper I concluded that in IR- luminous 
mergers the near-IR observations are dominated by the 
presence of a luminous, rotating young central stellar 
disk (YCSD) which contains a population of Red Su- 
pergiant (RSG) or Asymptotic Giant Branch (AGB) 
stars. Stellar disks have been directly detected in Hub- 
ble Space Telescope (HST) observa tions of the archetypal 
merger NGC 725 2 (iWhitmore et al. 1993) and the LIRG 
merger NGC 34 (|Schwe izer fc Seitzerii2007l ). with diam- 
eters ranging from 2-6 kpc. These are consistent with 
the extent of the excess light at K and //-band (RJ04, 
iHaan et al.!l2011[ ) ; the ob served size of rotati i ig molecular 
gas d i sks in mergers fe.g.FPupraz et al.!ll99dl : lWang et al.! 
II99I! : iDown es fc Solomon* '1 998[) and nume rical simula- 
tions (e.g. pBarncs 2002; Hop kins et aIll2008D . The rotat- 
ing YCSD affects the CTo.co measured in the centers of 
IR-luminous galaxies, which in turn affects the derived 
values of Moyn- However, at /-band, the presence of dust, 
which is more centrally concentrated due to the star- 
burst, behaves like a coronagraph. It masks the bright 
YCSD so that (Xo^caT refiects only the random motions 
of the old stellar population, probing the galaxy's true 
Afoyn- Figures 13-15 in Paper I demonstrated that the 
red (/ — K) colors within the central 1.53 kpc could be 
best explained by many magnitudes of extinction. While 
at larger radii the (I — K) colors (Figure 13 and 14 in Pa- 
per I) were consistent with the average colors observed in 
elliptical galaxies, supporting the proposition that dust 
is centrally concentrated in IR-luminous mergers. Thus, 
when viewed at near-IR wavelengths LIRG mergers ap- 

* The dividi ng line between radio loud and quiet at 6 cm is 10^'^ 
WHz-l Sr-l UMiller et al . 111990 : IStocke et al. 111993 : 1 Hooper et al.l 
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pear to have young stellar populations with Muyn < "m* , 
wher e m* is the stellar (not dynamical) m ass ~ 3x10^° 
Mq (|Blanton et alJ[200a IBell et alJ[2003l) . While at I- 
band they appear to have older stellar populations and 
Moyn > m*. 

The results in Paper I provide strong motivation for re- 
visiting the S88 paradigm of whether ULIRGs are mas- 
sive enough to form gEs and host QSOs. This paper 
presents the first results for 8 ULIRGs (part of a larger 
survey) using CTo from the CaT stellar absorption lines in 
conjunction with rest-frame /-band imaging. We probe 
two important questions: 1) Does the cr-Discrepancy ex- 
tend to the more luminous ULIRG population? and 
2) At /-band are the dynamical properties of advanced 
ULIRGs consistent with gEs, including the host galaxies 
of QSOs? 

All data and calculations in this paper assume Ho = 
75 km s~^ Mpc~^ and a cosmology of V,m — 0.3, V,\ = 
0.7 (qo = -0.55). Ah photometric results are in VEGA 
magnitudes. In this work, ULIRGs are strictly defined as 
LiR > 1012 °/q. LIRGs are strictly defined as 10" °iQ 
< iiR < 10"-^9Lq. 

2. SAMPLES 

2.1. ULIRG Sample 

The 8 ULIRGs analyzed in this paper were randomly 
selected (based on observability and available rest-frame 
/-band imaging only) from a larger, complete, volume 
limited (z < 0.15) sample of 40 advanced objects taken 
from the Infrared Astronomical Satellite (IRAS) 1 Jy 
Survey (Kim & Sanders 1998,) and the IR AS Revised 
Bright Galaxy Sample (ISanders et al.l[2Q03D . The IRAS 
1 Jy survey is a complete sample of 118 ULIRGs down to 
flux levels of /i/ = 1 Jy with Galactic latitude || &|| > 30°, 
declination S > -40° and 0.02 < z < 0.27. The Revised 
Bright Galaxy Sample is a flux-limited survey of galax- 
ies with a 60 /im flux density > 5.24 Jy covering the 
entire sky surveyed by IRAS. Late-stage ULIRGs were 
selected because CTo is u nlikely to ch ange substantially 
once the nuclei coalesce ()Mihosl[l999D . Based on numer- 
ical simulations and observations it marks the point at 
which the merger should exhibit properties in common 
with elliptical galaxies. HST F160W mCM0S2 images 
were used to confirm the presence of a single nucleus in 
each system (within the resolution limits of 59-182 pc). 
Six of the eight ULIRGs were observed with HST us- 
ing either the Advanced Camera for Surveys Wide Field 
Camera (ACS/WFC) or the Wide Field Planetary Cam- 
era 2 (WFPC2). Photometric data for the remaining 
two were obtained from the Seventh Release of the Sloan 
Digital Sky Survey (|Yorkll2000l : lAbazaiianll2009l hereafter 
SDSS DR7). Optical images of the ULIRGs are shown 
in Figure 1 of Appendix C. Table 1 lists the basic in- 
formation: Names, Right Ascension (R.A.), Declination 
(Dec), redshift (z). Log Lir,, and Galactic Reddening 
E(B - V). Lm is defined as the total flux from 8-1000 
(ISanders fc MirabelllT996l) using the four //Z^S* pass- 
bands (12, 25, 60 and 100 /xm). However, supplemental 
12 or 22/im photometry from the Wide-Field Infrared 
Survey Explorer was used in several cases where IRAS 
did not detect the ULIRG (see notes in Table 1). 

2.2. Comparison Samples 



2.2.1. SDSS i-hand DR7 Elliptical Sample 

In order to assess the signiflcance of the optically mea- 
sured values of cro,CaT and masses of ULIRGs, a compar- 
ison sample of elliptical galaxies was assembled from the 
SDSS DR7 which offers larger spectroscopic and photo- 
metric coverage over earlier releases and improvements 
in photometric and spectroscopic measurements. The 
comparison sample was selected to be volume limited (z 
< 0.15) to match the ULIRG sample. A total of 9,255 
elliptical galaxies were extracted from the SDSS using 
the selection criteria detailed in Appendix A. The se- 
lection criteria required that the elliptical galaxy must 
be present in both the photometric and spectroscopic 
databases. 

2.2.2. Radio Loud and Radio Quiet QSOs 

A relatively nearby (0.08 < z < 0.46) compari- 
son sample of 28 RLQ and 25 RQQ host galaxies 
was compiled from available photometry obtained with 
WFPC2 on HST and ground-based kinematic data. 
Only QSOs with confirmed ellipt ical host morpholog y 
were selected from th e samples of iBahcall et al.l (fl997l): 
iHamihon et al.' C2002'):'Dunl op et al.l ()2003D : iFlovd et al.l 
j2004); Hamihon ct al. (20Q1). The source papers ah 
note that their samples were designed so that: 1) the 
RLQ and RQQ subsamples are matched in terms of op- 
tical luminosity; 2) My < -23, representing L> L* galax- 
ies and ensuring that QSOs were selected; and 3) at z ~ 
0.4 the resolution of the WFPC2 cameras were sufficient 
to separate host from nucleus. In all cases the authors of 
the source papers performed extensive point spread func- 
tion (PSF) modeling using separate stellar observations 
to properly subtract the nucleus from the host. Although 
the assembled comparison sample of QSO host galaxies 
is heterogeneous, each of the source papers have demon- 
strated that their samples are statistically representative 
of the local QSO host galaxy populations. Moreover, 
the selection criteria employed by the source papers are 
rem arkably simflar wi th significant overlap. The D03 
and lFlovd et all (I2004D (hereafter F 04) samples are sub- 
samples from iMcLure et al.l (|1999D at 0.1 < z < 0.35, 
while the sample from Hamilton includes nearly all of 
the D03 sources, selected so that 0.06 < z < 0.4. The to- 
tal integration times of all of the QSO observations were 
checked with the latest version of the WFPC2 exposure 
time calculator to ensure they were sufficiently deep to 
properly sample the underlying host galaxy. The ba- 
sic information for the comparison sample of RLQs and 
RQQs are listed in Columns 1-6 of Table Bl in Appendix 
B. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. Optical F814 W Images 

The optical F8I4W Rlter was selected because the same 
filter was used for the /-band study in Paper I. That 
paper demonstrated that the F8I4 W filter is the best 
compromise between observing the old stellar popula- 
tion, which is used to probe the total MDyn, and avoid- 
ing light produced by RSG and AGB stars. Here, and in 
Paper I the F8I4W fAtei is simply referred to as /-band. 
The the mean differences among the ACS F8I4W, the 
WFPC2 F8I4 W, and the Cousins /-band fihers are less 
than a few hundredths of a magnitude. F8I4W images of 
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6 ULIRGs were obtained from the public HST archives 
(see Table 2 for more information). 

Five of the ULIRGs were observed with ACS/WFC 
as p art of the Great O bservatories All-Sky LIRG Sur- 
vey (Arm us et all [20091) . ACS/WFC is comprised of 
two 4096x2048 pixel CCDs, each with a platescale of 
0".049 pixel-\ providing a field of view (FOV) - 202"x 
202". This FOV is large enough to observe each of the 5 
ULIRGs completely. There is a gap of 50 pixels (2". 45) 
between the two CCDs. The observations employed a 
two position dither to fill the chip gap rather than the 
more common CR-SPLIT (two images taken at the same 
position) which is better for cosmic ray (CR) removal 
but leaves the gap with no data. As in Paper I, the 
ACS/WFC data were processed manually using the Im- 
age Reduction and Analysis Facility (IRAF) and The 
Space Telescope Science Data Analysis System (STS- 
DAS), which is a software package designed specifically 
for the reduction and analysis of HST data that works 
with IRAF. Individual exposures which have been cal- 
ibrated and flat field corrected were obtained from the 
archives for the 5 ULIRGs observed with ACS/WFC. For 
each ULIRG, the STSDAS task MULTIDRIZZLE was used 
to assemble individual dithered frames into a final mosaic 
image corrected for: geometric distortions and CRs; bad 
pixels set to a value of zero; and rotated to a position 
angle (P. A.) of 0°. This differs from the final drizzled 
image produced by the archive pipeline. Three bad pix- 
els were found in the center of IRAS 05189-2524 due to 
CR hits and warm pixels, and not due to saturation from 
over-exposure. The IRAF pixel editing task IMEDIT was 
used to replace the zeroed pixels with the values inter- 
polated from the surrounding pixels. The central region 
of IRAS 12540-1-5708 was found to be saturated. Both 
the diffraction spike generated from the bright core and a 
12x14 pixel rectangle in the center of IRAS 12540+5708 
were saturated. These were flagged by MULTIDRIZZLE 
and set to a value of zero. No images with shorter ex- 
posures were available to replace the flagged pixels and 
they were ignored in the subsequent analysis. Due to 
the two position dithering scheme used, MULTIDRIZZLE 
was unable to flag and remove a large number of CRs, 
particularly within the chip gap. Because the targets 
were centered in the ACS/WFC FOV, the chip gap runs 
through or close to the outer regions of the ULIRGs. As 
a result an algorithm was developed to remove CRs and 
is detailed in Appendix D. 

One ULIRG, IRAS F02021-2103, was observed with 
WFPC2, which is comprised of four 800x800 pixel 
CCDs. Three of the chips (WF2, WF3, and WF4) have 
a platescale of 0".099 and the fourth (PC) 0".046. This 
creates a non-symmetric FOV with a gap in coverage in 
the upper right quadrant. The observations for IRAS 
F02021-2103 were centered on the WF3 chip. The ob- 
servations were reduced using calibrated and flat field 
corrected WFPC2 science images obtained from the HST 
archives and processed with the STSDAS tasks WARMPIX, 
which fixes hot pixels, and CRREJ, which removes CRs 
and combines multiple frames into a single image. Geo- 

® IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



metric distortions were corrected by mu ltiplying the CR 
cleaned image with a correction image (iHoltzman et al.l 
[19951) . Finally, the image was trimmed to remove pixels 
vignetted by the pyramid shaped beam-splitter mirror. 

3.2. Near- Infrared F160W Im ages 

The F160W filter was selected because no similarly 
deep if-band data were available (as used in Paper I) 
and because the CO band-heads at 1.6 /^m were used for 
many objects to measure cto.co- Published or archival 
data were used only for the 5 ULIRGs with published 
values of (To,co- All observations used NICMOS with 
the NIC2 camera, a 256x256 pixel HgCdTe array with 
0".075 pixel-i platescale (19". 2 FOV). Photometric data 
for 3/5 ULIRGs were obtained from the literature. IRAS 
17208-0014 and IRAS 23365+3604 were analyzed from 
HST archival data. 

The NIC2 data were processed manually using IRAF 
and STSDAS. The raw (data received directly from the 
spacecraft without processing) files were used rather than 
the archive processed data in order to properly account 
for: 1) the presence of the NIC2 coronagraph; 2) the 
presence of bias jumps between the quadrants; and 3) 
the presence of electronic "bars" which appear as vertical 
stripes. The coronagraph shifts position over time and 
the anomalies vary with time requiring the individual 
raw frames to be processed manually. First, the STS- 
DAS task CALNICA was used to subtract dark current, 
correct for detector non-linearity, flatfield, convert to 
count rates, and identify and reject CRs. Next, the STS- 
DAS task PEDSUB was used to correct for bias jumps be- 
tween quadrants. This differs from the standard archival 
pipeline reduction which uses PEDSKY to remove both bias 
jumps and sky background. Because the objects fill most 
of the NIC2 array, using PEDSKY will result in a non- 
uniform over-subtraction of the background. In cases 
where electronic bars and other anomalies were found in 
individual raw frames, the STSDAS task NICPIPE was 
used instead. It allows the user to apply some or all 
of the steps from CALNICA. In this case, all steps except 
flat-fielding, conversion to count-rates, and CR rejection 
were applied to the data. Next the data were processed 
with BIASEQ, which corrects for drift in the bias levels 
during the course of MultiAccum exposures. The data 
were then processed through NICPIPE again, this time 
applying flatfields, conversion to count-rates and cosmic- 
ray rejection. The position of the coronagraph was de- 
termined in each individual frame and masked using the 
IRAF task IMEDIT. Other bad columns and hot pixels 
not removed with CALNICA or NICPIPE were manually 
identified and masked with the IRAF task IMEDIT. The 
frames were then processed with MULTIDRIZZLE in the 
same manner as the ACS/WFC F81 4W daXa, above, pro- 
ducing geometrically corrected images rotated to a P.A. 
of 0°. A comparison between this method and data pre- 
processed through the standard archive pipeline showed 
a significant improvement in signal to noise (S/N), in- 
cluding the detection of faint tidal features which would 
otherwise not be visible. 

3.3. Spectroscopy 

The optical spectra for all of the ULIRGs presented 
here were obtained with the Echellete Spectrograph and 
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Imager (jSheinis et al.|[200^ ESI) in echelle mode at the 
W.M. Keck II 10m observatory. Echelle mode employs 
a 20" long slit and cross-dispersed spectra with simul- 
taneous coverage of 0. 3927-1. 1068/im projected onto a 
2048x4096 pixel CCD. ESI has a fixed spectral resolu- 
tion of 11.5 km s^^ pixel^^. The final spectral resolu- 
tion scales with slit width. A 1".0 slit width (6.49 pix- 
els) was used for 7/8 ULIRGs. This corresponds to R 
~ 4000 or - 75 km s"^ A 0".5 slit width (3.24 pix- 
els) was used for IRAS 12540-1-5708. This corresponds 
to i? ~ 8300 or ~ 37 km s~^. In this paper only spectral 
orders containing the CaT stellar absorption lines (or- 
der 6 or 7 depending on redshift) were used. The scale 
along the spatial axis for order 6 and 7 are 0".168 and 
0".163, respectively. The integration time and P. A. for 
each ULIRG is listed in Table 2. Calibrations, including 
internal fiats and Hg-Ne, Xe, and CuAr arcs were taken 
at the beginning and end of the night. No changes were 
detected between fiats and arcs taken at the start and 
end of night. ESI spectra for three ULIRGs (IRAS 05189- 
2524, IRAS F10 378-H108, and IRA S 11387-^4116) were 
first presented in iRupke et al.l (|2002l ) but did not include 
CaT fj measurements. 

The data were reduced with IRAF. The reduction of 
the data and spectral extraction method used are nearly 
the same as those described in Section 3.1 of Paper I (e.g. 
correction to heliocentric rest velocities, spectrophoto- 
metric correction to remove instrumental signature and 
provide approximate fiux calibration, continuum normal- 
ization, and generation of the error spectrum) . The few 
differences with the methodology used in Paper I are de- 
scribed here. Due to the redshifts of the ULIRGs the 
CaT lines lie at wavelengths within a spectral region 
dominated by strong night-sky emission lines. This was 
corrected using the IRAF task BACKGROUND which mea- 
sures a sky spectrum at both edges of the slit and fits 
it with a polynomial (in this case a 1st or 2nd order 
Chebyshev polynomial) which is then subtracted from 
the spectrum. The polynomial fitting and subtraction 
is carried out column by column. In one case, the red- 
shift of IRAS F02021-2103 placed the CaT absorption 
lines coincident with some telUuric absorption lines. The 
IRAF task TELLURIC was used to correct for the presence 
of the telluric absorption lines using a featureless spec- 
trophotometric standard star. 

The spectra of 5/8 ULIRGs were extracted in a met- 
ric aperture of diameter 1.53 kpc (see Table 3). This size 
was selected to remain c onsistent with the l it erature (e.g. 
Jorgensen et al.l '1 9951 iSmith et al.l 11995': ' Pahre et al] 
1998a; Pahrc 1999'; iRothberg fc Jos eph 2006iB, herT 
after RJ06b, and Paper I.) However, in order to max- 
imize S/N, larger diameter apertures of size 3.77 kpc 
and 4.05 kpc were used to extract the spectra of IRAS 
F10378-I-1108 and IRAS 11387-^4116, respectively. Due 
to the brightness of the IRAS 12540-^5708 nucleus, the 
spectrum was extracted in a 0.82 kpc diameter aper- 
ture 2.08 kpc northwest of the nucleus. This is sev- 
eral times larger than the seeing of the observations, 
well within the reg ion dominated by stellar continuum 
(iDavies et al.l I2004D and the same method used by 
iTacconi et ail (|2002D : iDasvra et ahl (12006^ to measure 
(To, CO- Finally, the IRAF task CONTINUUM was used 
to identify and remove residual sky lines present from 
imperfect background subtraction in IRAS 17208-0014, 



IRAS 19542+1110, and IRAS 23365-1-3604. 

4. DATA ANALYSIS 

4.1. Photometry 

4.1.1. Galactic Reddening 

Photometry for all ULIRGs {HST and SDSS), com- 
parison ellipticals and QSO host galaxies have been cor- 
rected for Galactic reddening using E{B — V) values from 
iSchlegel et all (I1998D as presented in NEE|3 and assume 
R\r = 3.1 ()Fitzpatricklll999( ). These values were then 
scaled to the appropriate photometric filters: F8I4W 
and z for the ULIRGs; and the native filters listed in 
Table Bl for the QSO host galaxies using values from 
(jSchlaflv fc Finkbeineil 120111) . Scaling factors for the 
F675W a.iid F791 filters were co mputed using the York 
Extinction Solver (|McCalll [200l . The values used are 
listed in Tables 3, 4, and Bl. 

4.1.2. Measured Global Photometric Parameters 

Photometry was performed on the F8I4W images from 
ACS/WFC and WFPC2/WF3 and the F 1 60 W images 
from NIC2 to measure the global photometric parame- 
ters: the effective (or half-light) radius i?cff measured in 
metric units of kpc, the mean surface brightness within 
the effective radius (</i>cff ), and the total absolute mag- 
nitude (Ml or Mpieow)- The fluxes were measured in 
circular isophotes with fixed centers using the STSDAS 
task ELLIPSE. The position of the galaxy centers in the 
F8I4W images were determined from F160W NIC2 im- 
ages. As in Paper I, there were several cases where the 
nucleus was clearly visible in the F160W images but 
obscured or partially obscured in the F8I4W images. 
Foreground stars, bad pixels, artifacts or CRs missed 
by MULTIDRIZZLE were masked with a bad pixel mask 
created using the technique described in Appendix D. 
Masked items were set to a value of zero and ignored in 
the isophote fitting and flux measurements. An r^/"* de 
Vaucouleurs profile was fit to the isophotes produced by 
ELLIPSE for each galaxy. These surface brightness pro- 
files are plotted in Figure 2 of Appendix C. The angular 
effective radius (in arcseconds) from the de Vaucouleurs 
profile was converted to i?cff using the angular diameter 
and co-moving distance for our preferred cosmology. The 
values of </x>efi were derived from the 1^^^ pr ofile fits 
and were corrected for cosmological dimming ([Tolmanl 
Il930f ). The total Mi or Mpieow were computed by ex- 
trapolating the best-fit de Vaucouleurs model beyond the 
measured data and using the luminosity distance for our 
preferred cosmology. 

Table 3 lists the Mi, Rcb and </ii>cff values for the 
ULIRGs in the F8I4 W filter. Table 4 lists the mea- 
sured values of RcS at F160W and Mpigow for IRAS 
17208-0014 and IRAS 23365+3604. F160W photomet- 
ric results for IRAS F02021-2103, IRAS 05189-2524, 
and IRAS 12540+57 08, in Table 4 were obtained from 
iVeilleux et all (|2006D . For each of these ULIRGs, the 
sources were checked to ensure that an r^/^ profile was 
the better fit to the surface brightness profile. The val- 
ues of F160W RcS were converted to equivalent radii, 

The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, Cahfornia Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 
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which is the equivalent circular radius for measurements 
originally made with elliptical isophotes. The equiva- 
lent radii were computed using the semi-major and semi- 
minor axes, or the semi-major ax is and ellipticity (e.g. 
iMilvang- Jensen fc Jorgense 

4.1.3. Global Photometric Parameters of ULIRGs and 
Elliptical Galaxies from the SDSS 

No data reduction was performed on the actual SDSS 
DR7 images. Photometric data for IRAS F10378-f 1108 
and IRAS 11387+4116 were extracted from the SDSS 
DR7. The Sloan z-band filter was selected over the Sloan 
i-band filter based on two factors: 

1) At the redshifts of IRAS F10378+1108 and IRAS 
11387-1-4116 the Sloan z filter overlaps the rest- frame 
wavelength coverage of the F814 W filter; 

2) Because the wavelength coverage of the Sloan i-filter is 
bluer than the Cousins / and F8I4W Siteis, at these red- 
shifts fiux from strong emission lines of [01], Ha, [Nil], 
and [SII] could contribute significantly to the total ob- 
served flux for each galaxy. At the same time, the ob- 
served wavelength coverage of the z-band filter is now 
blue- ward of possible contamination from [SIII] at 0.9069 
and 0.9531 fxm. 

The metric equivalent radius for Res was computed 
from the SDSS parameters devRad, the half-light semi- 
major radius measured from a de Vaucouleurs fit to the 
galaxy light; and devAB, the axis ratio from the de Vau- 
couleurs best fit profile. Criteria 7 in Appendix A, that a 
galaxy is better represented by a de Vaucouleurs profile 
rather than an exponential profile, was verified for the 
two ULIRGs. Values of <^z>cff were computed from 
the SDSS parameters devMag, the total apparent magni- 
tude measured from the de Vaucouleurs fit to the galaxy 
light convolved with a double- Gaussian fit to the PSF; 
and the angular equivale nt radius Tcb (e.g. equation 7 in 
iHvde fc Bernardill2009bl hereafter HB09) and includes 
corrections for cosmological dimming. was also com- 
puted from devMag using the luminosity distance for our 
preferred cosmology. 

The devMag fiuxes were converted from AB magnitudes 
to VEGA magnitudes using the task CALCPHDT which is 
a part of the SYN PHOT (synthetic photometry) pro- 
gram in STSDAS (IHornd nm§. IKoornneef et~aLl [19861: 
iLaidler et "alll2005[) . CALCPHOT calculates synthetic pho- 
tometry for any input spectra or blackbody curves us- 
ing any filter transmission curve. It can be used to cal- 
culate fc-corrections and transformations between filters. 
A transformation value of (zvega — zab) = -0.51 was 
used. Due to large variations in the rest-frame ultraviolet 
and optical spectral energy distributions (SEDs) the k- 
corrections for ULIRGs at z ^ 0.1 can vary by A ifc 0.01 - 
0.09 mags fe.g. lSurace et al.lll998l : rTrentham et"alll 19991 ). 
Therefore, no ^-corrections were applied to the two 
ULIRGs. Although the same spectral regions are covered 
by the redshifted SDSS z filter and the rest-frame F8I4 W 
filter, the shape of the filter transmission curves for the 
two filters are very different. To assess the impact of this, 
we tested CALCPHDT with 12 sp ectral templates: a n el- 
liptical, SO, Sa, and Sb galaxies (jKinnev et al.lll996l ): six 
different starb urst galaxies with va riations in the amount 
of extinction ([Calzetti et al.|[r994[ ). including extinction 
similar to those observed in ULIRGs; and a composite 
spectrum of the Seyfert 2 NGC 1068, which includes ul- 



traviolet and optical lin es, nebular and po wer-law con- 
tinuum and cool stars (jLaidler et al.ll2005[ ). CALCPHOT 
yielded {F8I4W - z) = 0.19 ± 0.01 and 0.18 ± 0.01 for 
IRAS F10378-I-1108 and IRAS 11387-1-4116, respectively 
The small dispersion among the different templates is due 
to the similarity of the spectral features within the rest- 
frame F8 14 W wavelength range. These additional trans- 
forms were also applied to the values of devMag. For the 
remainder of the paper the observed Sloan z values for 
IRAS F10378+1108 and IRAS 11387+4116 will be re- 
ferred to as /-band. Table 3 lists the computed Mi, Rcs 
and </ii>cff values for IRAS F10378+1108 and IRAS 
11387+4116. 

Global photometric parameters for the SDSS ellipti- 
cals were computed in a similar fashion using the SDSS 
DR7 i-band values for devMag, devRad, devAB, red- 
shift, and extinction values, fc-corrections and a filter 
transformation from SDSS i-band to F8I4 W filter were 
computed using CALCPHOT and an e lliptical galaxy tem- 
plate spectrum ([Kinnev et al.l [19961) . The fc-corrections 
ranged from 0.03-0.16 mags with a filter transform of 
(FiSi^ 1/FvEGA ~ «abmag) — -0.59. Due to the large size 
of the SDSS elliptical sample, the computed photometry 
is not presented in the paper, but the selection criteria is 
provided in Appendix A and can be used to retrieve the 
sample from the SDSS DR7. 

4.1.4. RLQ & RQQ Host Galaxies 

The details of the data reduct fon methods used fo r 
these galaxies can be foun d in iBahcall et al.l (Il997f): 
[Hamilton et~al[ ([200 2): Dunl op et al.l (|2003[ ): lFlovd et al.l 
(|2004[ ): [Hamilton et al, t^OM)- The samples published 
in these papers were observed with WFPC2 on HST us- 
ing the F606W, F675W, F702W, F791W, or F8I4W fl- 
iers. The source papers provide surface brightnesses at 
the effective radius (/icfi), not surface brightness within 
the effective radius (</Lti>eff), and absolute and appar- 
ent magnitudes of each host galaxy and nucleus or PSF. 
In the source papers, these values were transformed from 
their observed filters to either rest-frame Johnsons V or 
Cousins /?-band magnitudes. To avoid adding additional 
uncertainties to the analysis, only the apparent magni- 
tudes of the host galaxies in the original HST filters were 
used. They were fc-corrected and transformed to rest- 
frame F8I4 W magnit udes using C ALCPHOT with an ellip- 
tical galaxy template ([Kinnev et al.l il996). An elliptical 
template was selected laecause all of the QSO host galax- 
ies have confirmed elliptical galaxy morphologies and ob- 
servations of the host galaxy spectra for nearly half of the 
sample indicate they are dominated by the presence of 
an older, quiescent stellar population at optical wave- 
lengths (Hu ghes et al. 2000; Canalizo & Stockton 2003 
[Nolan et al.l[2001l : ILetaw e et al. 2007; Wold ct al.ll2010i r. 

With the exception of .Bahcall et al. (1997.) the pub- 
lished values of ReS for the QSO host galaxies were mea- 
sured from elliptical isophotes. These were transformed 
to equivalent radii using ellipticities or semi-major and 
minor axes from the source papers and Hamilton {pri- 
vate communication) and converted to metric values of 
Rcs using our preferred cosmology. These, along with 
the transformed F8I4 W apparent magnitudes were used 
to compute <^i>cfi in the same way as for the SDSS 
galaxies. Mj, R^fi and <^i>cq are listed in Columns 
8-10 of Table Bl. 
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4.2. Velocity Dispersions 

4.2.1. Measurement of a o 

The details of the method used to measure cTq from 
the extracted one-dimensional ESI spectra are given in 
RJ06a, RJ06b and Paper I. The IDL routine VELOCDISP 
described in those papers and used here is based on 
a pixel-space direct fitting method to measure tio- 
This method is sim ilar to the technique described in 
iRix k Whitl (|1992| ). Briefly, the template stars are con- 
volved with a Gau ss-Hermite Polynomial, whic h is a 
modified Gaussian (Ivan der Marel fc FranxjUOQl . The 
18 template stars used for the fitting range from Gl to 
M7.5 giants to Kl to M5 supergiants (see Table 4 in Pa- 
per I) and are the same stars used in Paper I for the CaT 
and 2.3 /xm CO spectra. The fitting function has five pa- 
rameters: the line strength (7), which measures the ratio 
of the equivalent width of the galaxy to that of the tem- 
plate star; the mean recession velocity (uo), the central 
velocity dispersion (tJo, defined as ti in a 1.53 kpc met- 
ric aperture), the skewness (/la), and kurtosis (/14). The 
last two parameters characterize the departures from a 
Gaussian shape. The parameters are simultaneously fit 
to the data over the wavelength range 0.8480-0.8690 fim 
(|Barth et al. 2002). Bad pixel masks were used to mask 
out strong emission lines or regions of imperfect back- 
ground subtraction. Table 3 shows the best-fit results 
for the derived ctq, heliocentric recessional velocity (Vq) 
and best-fit template star for each ULIRG. The errors 
in Tables 4 were calculated by VELOCDISP using the er- 
ror spectrum for each galaxy. A more detailed discussion 
of this method of error analysis and a comparison with 
Monte Carlo simulations can be found in RJ06a, RJ06b, 
and Paper I. Figure 3 in Appendix C shows the CaT 
spectra for the eight ULIRGs (solid line), over-plotted 
with the best-fit convolved stellar template (gray dashed 
line). 

4.2.2. Aperture Corrections 

In order to reduce possible errors introduced by mea- 
suring kinematic properties in different aperture sizes, 
the values of CTo reported in this paper are either initially 
measured in a common aperture diameter of 1.53 kpc or 
corrected to this aperture size (see Tables 3, 4, and B2). 
This also applies to all published velocity dispersions, 
including the QSO host galaxies, SDSS ellipticals, and 
the CO velocity dispersions of ULIRGs. The c o rrecti ons 
used the scaHng relation from IJorgensen et al.l (|1995D : 



log 



(7 (do) 



a log — 

do 



(1) 



where a — -0.04^. do = 1.53 kpc and d is the metric 
diameter of the circular aperture. This is directly appli- 
cable to the SDSS ellipticals measured in a 3" diameter 
circular fiber. In the case of slit spectroscopy for the 
QSOs and ULIRGs, d is computed via: 

d ~ 1.025 X 2 X n (2) 

where x and y are slitwidth and extraction aperture, n 
is the metric scale (pc or kpc) in 1" for the galaxy com- 
puted from the angular diameter. The weighted average 

RFIO has a typographical error that reads a = 0.04 instead 
of -0.04 



radius along the summed portion of the slit for e ach QSO 
was taken from Table 3 in lWolf fc SheinisI ()2008D and the 
slitwidths for the CO observations of the ULIRGs from 
obtained from the relevant source papers (see Table 4). 

5. RESULTS 

The results here first address whether the a- 
Discrepancy extends to ULIRG luminosities, as sug- 
gested by Paper I. We then test the predictions of the 
S88 paradigm by using a combination of rest-frame I- 
band photometric data and (To,CoT to assess whether the 
global dynamical properties of the ULIRGs are consis- 
tent with those of massive ellipticals, including the host 
galaxies of QSOs, and whether they exhibit the same 
significant differences between optical and near-IR prop- 
erties. 

5.1. Extending the a -Discrepancy to ULIRGs 

A key result from Paper I was a demonstration that 
in advanced LIRG mergers, tTo,caT is systematically 
larger than (To,co- The same systematic discrepancy 
w as not observed in bon afide elliptical galaxies (Paper 
I. Vanderbeke et al.ll2011|) . Paper I further posited that 
the CT-Discrcpancy should also extend to ULIRG lumi- 
nosities. We now test this hypothesis by comparing the 
values of CTo.CaT with published values of (To,co- Figure 
1 shows the five ULIRGs in the sample which have pub- 
lished values of Uo.co (see Table 4) compared with their 
""o.CaT presented in this paper (see Table 3). The dotted 
line represents (To,caT = Co.co- Taking into account the 
errors, the fTo,CaT values plotted for the 5 ULIRGs lie 2.5- 
22.9(7 away from the expected (To.co values. The relative 
difference between (To.CaT and <To,co for each galaxy was 
characterized in Paper I by the parameter o-prac: 



CFrac 



Co, CaT — fo,CO 
Co,CaT 



(3) 



This parameter was then compared with Log Lir (Figure 
4 in Paper I) for both advanced mergers and bonafide el- 
liptical galaxies. Lir seems a natural starting point for 
comparison, given that the degree of the a-Discrepancy 
appears to be greater for LIRGs than non-LIRGs, and 
non-existent in the control sample of ellipticals (Figure 1 
of Paper I). The comparison between trprac and Log Lir 
was quantified using the Pearson Correlation coefficient 
(r) which tests the degree of linear correlation between 
two independent data sets. The value of r ranges from 
-1 to -1-1 (anti-correlation to perfect positive correlation). 
In Paper I, a comparison was made for the entire sample 
of advanced merger remnants (9.51 < Log Ljr < 11.96), 
yielding a value of r = 0.77 (a strong correlation). The 
comparison sample of ellipticals in Paper I showed no 
correlation (r — 0.06). In order to determine whether 
this correlation extends to ULIRG luminosities plotted in 
Figure 2 are aprac and Log Lir values for the 5 ULIRGs 
from this paper and all advanced mergers from Paper I. 
The correlation in Figure 2 is r = 0.75^° °^. The errors 
on this correlation were compu ted using a "jackknife" re- 
sampling method (|Tukevlll958l ) in which the computation 
of the Pearson Correlation for the sample of n pairs of 
data points is done using n - \ pairs of data points. This 
allows for n computations of the Pearson Correlation to 
be made and a standard deviation to be computed to 
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test the robustness of the correlation (i.e. that one point 
may drive an apparent correlation). 

The algorithm FITEXY was applied to the galaxies in 
Figure 2. It employs a double-weighted least-squares 
(DWLSQ) fit using a minimization method that ac- 
counts for errors in b oth variables (i Press et al.i il992i : 
iFeigelson fc Babulll992[ ). The result is: 

CTFrac = 0.18±° "i Log iiR - 1.79±" " (Log LiR > 9.5) 

^ (4) 

This is similar to Equation 7 in Paper Fn . The Log Lir- 
f Frac correlation first shown in Paper I extends to ULIRG 
luminosities, i.e. the relative difference between Co^caT 
and (To, CO grows with larger Ljr. 

5.2. The Dynamical Properties of ULIRGs 
5.2.1. The Fundamental Plane 

Figure 3 shows the /-band Fundamental Plane edge 
on (diagonal solid solid line). It is a two-dimensional 
plane embedded within the three- dimensional parame- 
ter sp ace of CTo, ReS, and <n>efi ()Diorgovski fc DavisI 
I1987D . Early- type galaxies lie on the FP, late- type galax- 
ies do not. A relationship similar to the FP can be de- 
rived fro m the V irial Theorem (e.g. 'Fa ber et al.l 119871 : 
|DjorgovskretaD[i988; Bender et al. 199^. The FP can 
be used as a diagnostic tool to probe whether the dynam- 
ical properties of a particular galaxy, or group of galaxies 
are similar to those of elliptical galaxies. In this section 
the two fundamental goals are: 1) to ascertain whether 
the advanced ULIRGs lie on or close to the FP; and 2) 
to compare where the ULIRGs lie relative to the QSO 
host galaxies. 

The FP plotted in Figure 3 is from the orthogonal fit 
in Table 2 of HB09 and was derived from ~ 50,000 early- 
type galaxies (z < 0.36) using photometry from SDSS 
Data Releas e 4 an d values of cr from SDSS Data Re- 
lease 6 (see iHvde fc Bcrnardi 2009a, for more details). 
As noted in Appendix A, we used the same parameters 
as HB09 to select the comparison sample of ellipticals 
but with tighter restrictions on morphology in order to 
select only ellipticals and exclude SO galaxies. Thus, it 
is a sub-sample of IIB09. 

The Sloan i-band was transformed to Ho — 75 and 
to the HST ACS/WFC F8I 4 W fil ter using the same 
method described in Section 14.1.31 In Figure 3, and 
subsequently for all comparisons, ULIRGs are plotted as 
filled circles, RQQs as open flattened diamonds, RLQs 
as open squares, and the SDSS ellipticals are plotted in 
four groups based on their values of CTo (see Figure 3 
caption). The first group roughly corresponds to sub-m* 
ellipticals, the second to ~ m* ellipticals and the last 
two to progressively more massive systems. This adds 
an additional dimension of information to these plots; 
and in the case of Figure 3 separates the two param- 
eters which comprise the x-axis, revealing a gradient of 
increasing CTq from left to right across the FP. This gradi- 
ent reflects a previously reported steepening of the slope 
of the FP as a functio n of CTo fe.g. iJorgensen et"aLlll996l : 
iBernardi et"ani2003al HB09). Also of note IS an appar- 

Equation 7 in Paper I was computed using a least-squares fit 
with errors in Y array. We have recomputed the fit using FITEXY, 
which produces crprac = 0.18±'' 0^ Log Lir - 1.78±°12 



ent thickening and slight warp or curvature of the FP 
at small i?eff and fow gp which is not due solely to ob- 
servational errors (IJorgensen et al.lll996t IBernardi et al.l 
l2003a; Nig oche-Netro et al.ll2009r iIB09V 

To determine if ULIRGs lie on or close to the FP, the 
scatter of the residuals or r.m.s. of the ULIRG sample 
relative to the FP is compared with those of the SDSS 
Ellipticals. By definition, these ellipticals lie on the FP 
because they are a sub-sample of the early-type galaxies 
used to construct the FP itself. The SDSS ellipticals have 
r.m.s. = 0.18 dex (in units of kpc), while the ULIRGs 
have r.m.s. = 0.29 dex. The residuals of the ULIRGs 
range from -0.72 to 0.30 (the SDSS ellipticals range from 
-0.68 to 0.40), with a mean and median of -0.05 and 
0.04, respectively, indicating no systematic offsets for the 
sample. In general, the ULIRGs lie well within the scat- 
ter of the SDSS ellipticals, with four lying on the FP 
(within their errors). We conclude that the ULIRGs lie 
on or close to the FP like bonafide ellipticals. We note 
that IRAS 19542-HlllO is 3.9ct from the FP, making it 
an outlier, although there are SDSS ellipticals which are 
similarly distant from the plane. 

Where do ULIRGs lie on the FP and how does their 
location compare with those of gEs, including QSO host 
galaxies? To quantify this, we used a variation of the 
Kolmogorov-Smirnov (KS) two-sided (i.e. comparison 
between two empirical distributions) te st which is ap- 
plicable to two-dimeiisional data sets (iPeacockl 119831 : 
iFasano fc FranceschinillT987l: IPress et al.l ll992y The KS 
test itself probes the Null Hypothesis that the two dis- 
tributions to be compared have the same distribution. 
It is a non-parametric test, meaning no assumption is 
made about the form o f the distribu tion except that it 
must be continuous (e.g.'Massey"195l!). A standard rejec- 
tion threshold of 95 % (also known as th e 0.05 confidence 
level) was selected (jFished IT9251 . 119901 ) for the analysis. 
If the Null Hypothesis can be rejected at a greater con- 
fidence it will be stated, otherwise a statement of rejec- 
tion or non-rejection will always refer to the 95% level. 
The two-dimensional form of the KS two-sided test was 
designed to test the empirical distribution of data points 
on a plane and provide a goodness-of-fit statistic without 
the problems which arise from binning (i.e. test) or 
assumption of a particular shape to the distribution. The 
te st statistics were co mputed using the methods outlined 
in IPress et al.l ()1992| 1 which arc in turn, based on modifi- 
cations to the KS statistic (jStephens 197Q). These allow 
for computation of the 2 D statistic beyond the limite d 
case examples provided in lFasano fc Franceschinil ()1987D . 
We note two important caveats for the KS t ests u sed 
here and throughout the remainder of Sections 15.2.11 and 
15.2.21 First, the reported errors for the RLQ values of 
(To listed in Table B2 are significantly larger (± 17-34%) 
than those of the RQQs, ULIRGs, or SDSS ellipticals. 
They are also significantly larger than those of any other 
parameters examined in this paper. Such large errors 
may affect the KS tests. A set of 10,000 Monte Carlo 
simulations were performed in which a new value of CTo 
was randomly generated from within the range of Uo ± 
A(To for each RLQ. The KS test was then re-run for each 
Monte Carlo simulation to check for changes in the re- 
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sultJ*^. As an aid to the reader, the resuhs of aU KS 
tests (ID and 2D) are summarized in Appendix C in Ta- 
ble CI. The values in parentheses in Table CI show the 
percentage for which the results reamin the same in the 
simulations compared to the actual test result. Second, 
although IRAS 19542+1110 is a 3.9cr outlier on the FP, 
excluding it from the KS tests does not change the re- 
sults presented here and throughout Sections 15.2.11 and 
[5X21 

The Null Hypothesis cannot be rejected when the 
ULIRGs are compared with: the 6 RLQs and all 8 of 
the QSO host galaxies (RLQs and RQQs). The Null Hy- 
pothesis can be rejected when the ULIRGs are compared 
with the entire distribution of the SDSS sample. When 
the comparison is restricted to the SDSS ellipticals in the 
two largest bins, 165-225 km s^^ and 225-420 km s"'^, the 
Null Hypothesis cannot be rejected for either bin. 

We now compare the velocity dispersions of the sam- 
ples and d efer comparisons for Log R^s and </j.i> to Sec- 
tion [5?2]4] where the RLQ and RQQ samples are signifi- 
cantly larger. A standard (ID) two-sided KS test com- 
parison was made for ctq between the ULIRGs and QSOs, 
and the ULIRGs and SDSS ellipticals. It should be noted 
that the methodology for the standard two-sided KS 
test uses the tables for small sampl es or equations for 
large samples originally published in iPearson fc HartlevI 
(|1972f ) as well as the modified KS test which is rou- 
tinely used in programming language libraries and sta- 
tistical software (e.g. Fortran, C, C-f+, IDL, Python, 
R) for comparisons among samples of any size with- 
out the need for comprehensive tables (jStephensI 119701 : 
iPress et aLlll992D . If the two methods disagree, it will be 
noted, otherwise it is assumed that both methods yield 
the same result. 

The Null Hypothesis can be rejected when the distri- 
butions of (To are compared for the ULIRGs and RLQs, 
but only 26% of the time. The Null Hypothesis cannot be 
rejected when the ULIRGs are compared with all 8 QSO 
host galaxies (6 RLQs -I- 2 RQQs). The Null hypothesis 
can be rejected (and at the 99% level) when comparing 
the distributions of do for the ULIRGs and the entire 
SDSS comparison sample. 

Figure 3 and the 2D KS tests show that at /-band 
ULIRGs lie on the FP in a region where M >> m* 
ellipticals are found. This is in contrast to pure near- 
IR studies which showed ULIRGs are systematically off- 
set from the FP in regions dominated by low to in- 
termediate mass ellipticals (M < m*) and therefore 
could not be the progenitors of QSO hos t galaxies (e.g. 
Shier fc FischeilllgMlGenzel et al.| [200lt iTacconi et at] 
20021: iDasvra et al.ll2006[ ). The comparison with QSO 
host galaxies is less clear. The do distribution may not 
be consistent with that of RLQs. The uncertainty arises 
from the large errors associated with values of Uo for the 
RLQs. However, one cannot rule out similarities between 
ULIRGs and the QSO host galaxy population as a whole. 

5.2.2. Dynamical Masses & Stellar Populations at I-band 

Paper I demonstrated that the observed dynamical 
properties of LIRGs are different at /-band and /T-band, 

The same Monte Carlo tests applied to the other parameters 
and applied to the errors in (To for the ULIRGs and SDSS ellipticals 
yielded no change in the KS ID and 2D test results 



an effect not seen in elliptical galaxies. It showed that for 
a given LIRG, the presence of a central, relatively young 
population of RSG and/or AGB stars dominates the K- 
band light. As a result, MDyn measured at /sT-band is sig- 
nificantly smaller than Mnyn measured at /-band. The 
apparent effective ages also typically younger at /T-band 
(see Figure 9 in Paper I). At /-band, this population is 
largely obscured due to dust, permitting the kinematics 
of the older, more global population to dominate the ob- 
servations. Figure 4 is similar to Figure 9 in Paper I. 
It shows Moyn vs M/L at /-band. The ULIRGs, RLQs, 
RQQs, and SDSS ellipticals are plotted in Figure 4 (same 
symbols as Figure 3). L represents the total luminosity 
and the masses shown are the total virial A/oyn of each 
galaxy: 



Mnyn = ^ -^-^ (5) 



(IProvedal 119581 iFishI 11964 


iRood et all 119721: 


Tonrv & DavisI 119811 iBinnev 


119821: Bacon et al.l 


19851: iRichstone & Tremaind 119861: iMathewsl 119881: 


Bender et all 119890 where 


is the effective radius 



from the de Vaucouleurs fit, G is the gravitational 
constant, and k =6 (which takes into account the 
variations in shape, size, and inclination of spheroids; 
the impact of rotation on Uo] and that a varies with 
radius). The values for MDyn are listed in Table 3 for 
the ULIRGs and Table B2 for the QSO host galaxies. 
The vertical dotted line represents an m* elliptical 
galaxy. Overlaid are two sets of models representing the 
evolution of M/L for a single stellar populati on (SSP). 
The pair of models on the left are from (jMarastonl 
120051 hereafter M05). The sohd vector is an SSP with 
solar metallicity and a Kroupa initial mass function 
(IMF). Changing from solar to either half or twice 
solar metallicity causes only a slight shift (< 0.1 dex) 
up or down in M/L, respectively. The vector shown 
in light grey parallel to the Kroupa vector is an SSP 
with solar metallicity and a Salp eter IMF. On the 
right are updated SSP models from (jBruzual fc CharlotI 
[200l the models are hereafter referred to as CB07). 
The solid vector is an SSP with solar metallicity and 
Chabrier IMF and the light grey vector parallel to it 
is an SSP with solar metallicity and Salpeter IMF. 
Using sub-solar (0.4 Zq) or more than solar (1.5 Zq) 
metallicity decreases or increases the M/L values by 
no more than 0.15 dex. The M05 and CB07 models 
generally agree with each other at /-band, although the 
latter shows more variation/degeneracy in M/L at i ~ 
1-1.2 Gyr, while M05 shows some variation/degeneracy 
at 0.2-0.4 Gyr. These are likely related to differences in 
their treatment of thermally pulsing as^^nptotic giant 
branch (TP- AGB) stars (e g. see M05, iMaraston et aU 
[200l iBruzuall [20€7l [20m ). Such differences have a 
more pronounced impact on the analysis of the stellar 
populations of (U)LIRGs at near-IR wav elengths (e.g. 
lRothberd[2009l : iRothberg fc FischeiilMOal and Paper I) 
The horizontal placement of both vectors are for display 
purposes only. 

Figure 4 demonstrates that at /-band, all of the 
ULIRGs have MDyn > With the exception of 

IRAS 19542-f 1110, the ULIRGs lie at the upper end of 
the mass distribution. As noted earlier, ULIRGs are 
known to contain massive quantities of H2 and have 
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prodigious star-formation rates. Var i ous methods for 
estim ating the SFR (e.g. iKennicuttI 119981: lYun et al.l 
1200 ID indicate the 8 ULIRGs plotted in Figure 4 have 
SFRs ~ 100-500 Mq yr-i. Yet, when the M/Li values 
are compared with the SSPs, the stellar populations 
appear to be old and evolved for 7/8 ULIRGs. A 2D KS 
test comparison of the distribution of the ULIRGs with 
the galaxies in the /-band M-oy^-M/ L plane indicates 
that the Null Hypothesis can be rejected when the 
ULIRGs are compared with the QSO host galaxies 
(either the RLQs alone or all 8 QSOs). However, the 
large RLQ Uo errors weaken this result significantly, 
especially when the entire QSO sample is considered (see 
Table CI for details). In addition, the Null Hypothesis 
can be rejected for the entire comparison sample of 
SDSS ellipticals, although the Null Hypothesis cannot 
be rejected when the ULIRGs are compared with the 
SDSS ellipticals in the 165-225 km s'^ and 225-420 km 
s~-^ bins. . 

It should be noted that Figure 4 presents a simplistic 
approximation by assuming a single burst population 
is representative of the entire galaxy, i.e. all stars are 
of the same age. Moreover, the M/Li models shown in 
Figure 4 are for stellar masses only, whereas the plotted 
data are dynamical masses. Although SSP models may 
be an adequate approximation for a typical quiescent 
elliptical galaxy, the models do not take into account 
the presence of non-stellar matter (e.g. gas, dust, etc), 
nor multiple stellar populations (such as a younger 
population unseen at /-band, as we posit for LIRGs 
and ULIRGs). Both the additional ISM mass and the 
effects of extinction of the /-band flux from dust will 
increase the plotted values of M/ L. Thus, the dynamical 
M/ L is almost always greater than M/ L derived from 
stellar population models, and in the case of younger 
galaxies, the inclusion of two populati ons brings the 
mode ls more in line wi t h the d ata (e.g. iGerhard et al.l 
[200lt ICappellari et~all [20061: Ide Jong fc Belli 120071) . 
However, the main point we make here is that at /-band 
the ULIRGs appear closer in mass and age to older 
gEs (with masses well above m*) in contrast to results 
obtained for the same ULIRGs at near-IR wavelengths. 
We do not attempt to derive absolute stellar population 
ages for the ULIRGs using a single bandpass, we simply 
point out that the comparisons in Figure 4 appear to 
indicate little or no evidence for the presence of young 
stars at /-band, even though they are clearly present at 
other wavelengths. 

Separating the two parameters in Figure 4, we now 
focus exclusively on MDyn using a standard two-sided 
KS test. First, for the QSOs, the Null Hypothesis 
can be rejected (and at the 99% level) when the Moyn 
distributions of the ULIRGs are compared with the 
RLQs. Monte Garlo simulations show little change, 
except at the 99% level. The Null Hypothesis cannot be 
rejected when the Mnyn distributions of the ULIRGs are 
compared with the entire sample of QSO hosts. Next, 
the Null Hypothesis can be rejected (and at the 99% 
level) when the ULIRGs are compared with the entire 
sample of SDSS ellipticals. The only qualification is that 
the Null Hypothesis cannot be rejected for the largest 
(Jo bin. What these results show is that at /-band the 
ULIRGs are consistent with gEs possibly including RQQ 
host galaxies. They may not be consistent with RLQ 



host galaxies, but some doubts remain. 

5.2.3. The Dynamical Masses & Stellar Populattons of 
ULIRGs: Optical vs. Near-IR 

As noted earlier, previous dynamical studies of 
ULIRGs carried out almost exclusively in the near-IR 
had concluded that for these systems Moyn < m* ■ Figure 
5 demonstrates how the cr-Discrepancy in Section r5.1l can 
account for the differences between the results in Sections 
I5.2.1l and l5.2.2l and earlier near-IR only studies. This fig- 
ure presents a relative comparison of where ULIRGs lie 
in the Mny-n-M/L plane at /-band (top) and at F160W 
(^ //-band, bottom). Only the ULIRGs which have kine- 
matic and photometric observations at both wavelengths 
are shown. The ULIRGs are represented by letters in 
Figure 5, each letter corresponding to a specific galaxy 
(see the caption for Figure 5 and Table 4). The two 
plots share the x-axis (MDyn) in order to enable a direct 
comparison of the mass computed at both wavelengths. 
The F160W photometry and kinematic data from the 
near-IR CO band-heads are presented in Tabic 4. As in 
Figure 4, SSP models from M05 are shown with Kroupa 
(solid vector) and Salpetcr (light grey) solar metallicity 
IMFs on the left and CB07 for a Chabrier (sohd) and 
Salpeter (fight grey) IMFs on the right. The SSP vec- 
tors plotted in the F 1 60 W 'panel were transformed from 
the original M05 Johnson //-band values to F160W us- 
ing (// — F160W^) colors computed by processing the 
grid of M05 SEDs at each age with SYNPHOT. At F160W, 
the discrepancy between the M05 and CB07 models are 
more pronounced than at /-band, particularly for the 
M05 Kroupa IMF and the CB07 Chabrier and Salpeter 
IMFs. The relative ages appear to be offset by ^ 0.5 
dex between the M05 and CB07. However, regardless of 
the specific age, ULIRGs clearly appear to be older and 
more massive at /-band than at F160W. As first noted 
in Paper I, this effect is not seen in elliptical galaxies be- 
cause they do not have a young stellar population that 
dominates observations in any wavelength. 

5.2.4. The Kormendy Relation 

The analysis presented above is limited by the num- 
ber of QSO host galaxies with Uo measurements, partic- 
ularly RQQs. The </i>off-Log /?cff plane (also known 
as the Kormendy Relation) is a photomet ric projec- 
tion of the FP for early-type systems (e.g. iKormendvl 
119771 1198^ . and, like the FP, is independent of galaxy 
environment (e.g.. IPahre et al.l ll998bl: iReda et al.l l2004i: 
iNigoche-Netro et al.ll2007[ ). Although the Kormendy Re- 
lation (KR) has significantly more scatter than the FP, it 
is often used as a "cost-effective" proxy because the stel- 
lar absorption line spectroscopy needed to measure tJo 
can be time-consuming and/or difficult to obtain. This 
is especially true for QSOs because observations must 
contend with the effects of the bright nucleus which can 
swamp the underlying host galaxy. Until recently, nearly 
all studies investigating the dynamical properties of QSO 
host galaxies and their relationship to ellipticals, spirals, 
and mergers have relied solely on photometric observa- 
tions. Here, the KR is used to increase the number of 
QSO host galaxies for comparison from 6 to 28 RLQs, 
and from 2 to 25 RQQs. The main goals of this section 
are: 1) to ascertain the positions of advanced ULIRGs 
with respect to the KR; and 2) to determine where the 
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ULIRGs lie relative to the elliptical QSO host galaxies. 

Figure 6a {left) shows the /-band KR (sohd hne) de- 
rived from the entire comparison sample of SDSS ellipti- 
cals using the DWLSQ fitting method described earlier. 
The derived /-band KR is: 

< >^ff = 17.52±°°i + 2.26±°°iLog Raff. (6) 

with r.m.s. — 0.39 dex (in units of mag arcsec"^) and a 
Pearson Correlation Coefficient of r = - 71. T his is sim- 
ilar to the fit from Nigochc-Netr o et al.l (j2008D for 8,664 
early type galaxies at SDSS i-band (z < 0.36). Their 
slope and intercept is 2.52^" °^ and 17.84*° °^, respec- 
tively (transformed to F8I4W, the intercept is 17.25), 
with r.m.s. = 0.42 dex and a Pearson Correlation Coef- 
ficient of r = 0.72. 

The ULIRGs, RLQs, RQQs, and SDSS ellipticals are 
plotted in the left panel of Figure 6a (same symbols as 
Figure 3). The diagonal dotted shows the loc u s of an L* 
galax y (|Blanton et all I200I iBeh et al.l [200l iHill et al.l 
l2010r ). Also shown are the 1,2, and 3a dispersions of the 
KR (dotted grey lines). As in Figure 3, the SDSS ellipti- 
cals follow a gradient in their distribution when grouped 
by (To . The SDSS ellipticals with the largest CTq lie almost 
exclusively above the KR and the ones with the smallest 
(Jo he almost exclusively below it. In addition to fitting 
a KR to the entire SDSS comparison sample, fits have 
been made for each of the SDSS sub-samples (binned by 
(Jo) as well as the ULIRGs and QSO host galaxies. The 
number of objects in each sample, the coefficients of the 
fit, r.m.s. and Pearson Correlation Coefficient are pro- 
vided for each sample in Table 5. The slope of the KR 
does not vary significantly among the first three (Jo bins, 
but does decrease significantly for the bin with the high- 
est (Jo. This is somewhat d ifferent than the results from 
INigoche-Netro et all (|2008f) . which found that the slope 
of the KR changes significantly as a function of luminos- 
ity when binned in 1 mag intervals and in intervals of 
increasing luminosity. 

The r.m.s. of the 8 ULIRGs is 1.01, within 2.5a of the 
KR. However, IRAS 19542+1110 is a significant outlier 
from the /-band KR ('^ 6a). Excluding it, the r.m.s. 
of the ULIRGs decreases to 0.52 1.3a from the KR). 
Once again, a 2D KS test was used to compare the dis- 
tribution of the ULIRGs with all of the SDSS ellipticals. 
Here, the results do change when IRAS 19542+1110 is 
included or excluded. When included, the Null Hypoth- 
esis cannot be rejected. When it is excluded, the Null 
Hypothesis can be rejected at the 95% level. However, 
when the ULIRGs are compared with the SDSS ellipti- 
cals in the 165-225 km s'^ and 225-420 km s^^ bins, 
then the Null Hypothesis cannot be rejected, whether 
or not IRAS 19542+1110 is excluded. This implies that 
the ULIRGs are always consistent with the most massive 
SDSS elhpticals, and the inclusion of IRAS 19542+1110 
also makes them consistent with a broader range of el- 
lipticals. 

When the ULIRGs are compared with the RLQs, 
RQQs, and the QSO host galaxies taken together as one 
sample, the 2D KS test shows that in all cases the Null 
Hypothesis cannot be rejected. Although the majority 
of ULIRGs and QSOs fall within 3a of the KR, they lie 
systematically above the relation. The computed slopes 
and intercepts are listed in Table 5. Figure 6b {right 



panel) shows the best-fit KRs for the ULIRGs, RLQs, 
and RQQs from Table 5. The fits are plotted as shaded 
regions which account for the la errors in both slope and 
intercept. 

The ULIRG and QSO fits have steeper slopes than the 
SDSS ellipticals (including the sub-samples binned by 
(Jo). The fits to the ULIRGs and QSOs are consistent 
with each other. Figure 6b shows the considerable over- 
lap in KR param eter space among UL IRGs, RLQs, and 
RQQs. Similarly, l^meuxeFaD ([200^ found at F160W 
that the slope of the KR fits to PA'F-subtracted PC QSOs 
and ULIRGs (in which the bulk of the star-formation 
should have been removed in the PiSf-subtraction ) , along 
with QSOs from D03 and [Hamilton et all (|2008[ ). trans- 
formed assuming {R — H) = 2.8, were indistinguishable 
from each other, but still significantly steeper than in- 
active ellipticals. The QSO results in Table 5 and Fig- 
ure 6b also match results from earlier attempts to place 
QSO host galaxies o n the KR at optical wavelengths (e.g. 
lO'Dowd et aL|[200l D03,F04, see bottom of Table 5). 

A standard two-sided KS comparison of each individ- 
ual parameter in the KR (</i>off and Log Raff) between 
the ULIRGs and the RLQs, RQQs, and the two QSO 
samples together produces the same result in all cases: 
the Null Hypothesis cannot be rejected. Given that the 
ULIRGs are now compared with much larger samples 
of RLQs and RQQs, Figure 6 presents stronger support 
than the FP or Mj^y^-M/ L plane for the assertion that 
ULIRGs are consistent with the host galaxies of QSOs, 
both as a single population, and when divided into RLQs 
and RQQs. 

Finally, since the RLQ and RQQ photometric samples 
are sufficiently large, one can compare them with each 
other using both the 2D KS test and the ID compari- 
son for each parameter in the KR. Here, the results are 
quite interesting because the Null Hypothesis cannot be 
rejected for the 2D case nor can it be rejected for </x>efi. 
However, for Log Reff the comparison is less clear. Using 
the modified KS test, the Null Hypothesis can be rejected 
at the 95% level, but using the standard KS formula for 
large samples it cannot be rejected. Such a difference be- 
tween the methods suggests one should err on the side of 
caution in making any strong statements about whether 
the RLQs and RQQs are significantly different in a sta- 
tistical sense. Both D03 and F04 reached a similar con- 
clusion, noting that the mean R^ff of RLQs and RQQs 
are the same within the la errors. These results raise 
doubt about the intrinsic differences between RLQs and 
RQQs. This reinforces the need for more kinematic ((Jo) 
measurements for RLQs and RQQs. 

6. DISCUSSION 

The goal of this paper has been to use the first results 
of the /-band dynamical survey of advanced ULIRGs to 
address two key questions: 1) Does the a-Discrepancy 
extend to the more luminous ULIRG population? and 
2) At /-band are the dynamical properties of advanced 
ULIRGs consistent with gEs, including the host galaxies 
of QSOs? Here, we brieffy discuss the implications for 
the results presented so far. 

6.1. The a -Discrepancy 

The (j-Discrepancy does appear to extend to ULIRGs 
such that (Jo,caT > Co, CO- Moreover, the results here 
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agree with what was posited in Paper I, namely the 
correlation between Log Lir, and the cr-Discrepancy, 
and in turn, the predicted range ol (To,caT values for 
ULIRGs. Just as with the LIRGs in Paper I, the ULIRGs 
can be described as Janus-hke. Like the Roman de- 
ity, they present two different faces depending on how 
they are viewed. At /-band the face of an old stel- 
lar population is observed, while at near-infrared wave- 
lengths the face of a young stellar population domi- 
nates (see Figure 19 in Paper I). Since in Paper I, and 
in this work, we are observing only single nuclei merg- 
ers, the luminosity evolution of the disks in these sys- 
tems may be nearly monotonic and decreasing. In this 
subclass of (U)LIRGs, the lower the Lir, the further 
along will be important processes su ch as feedback that 
clears out the star-forming ISM (e.g. iFischer et al.l[201Cll : 
iFerugho et all 120101 : iSturm et al] |201lD and t"he subse- 
quent aging of the starburst. As the YCSD becomes 
fainter, its dominance in both photometric and kine- 
matic measurements in the near-IR subsides, reducing 
the observed cr-Discrepancy. As the lack of this dis- 
crepancy in bonafidc ellipticals demonstrates (Paper I, 
FVanderbeke et al. 2011) at some point in the evolution- 
ary sequence ao,opUcai = <To,near-iR- Part of this expla- 
nation is an oversimplification because it assumes that 
every merger reaches a ULIRG stage and that LIRGs 
are stages before and after the luminosity peak. On the 
other hand, regardless of whether all LIRG mergers will 
be or at some point have been ULIRGs, lower Ljr means 
a less luminous YCSD, and a less dusty nuclear region. 

6.2. The Evolution of ULIRGs into QSOs 

The extension of the a-Discrepancy to ULIRGs leads 
directly to the second question; whether the structure 
and kinematics of ULIRGs are consistent with those 
of gEs, including QSO host galaxies. We address this 
in two steps, beginning with a comparison to elliptical 
galaxies in general. Previous results from near-IR stellar 



kinematics and photometr y (e.g. Genzel et al.' '200^ 
Tacconi et al. 2002; Davi es et al.l 120041 : iColina et al] 
2005t iDasvra et al.l (20060 concluded ULIRGs were the 
progenitors of low- intermediate mass ellipticals (< m*) 
based on two arguments. The first was simply that the 
measured values of CTo obtained from near-IR stellar 
lines were significantly smaller than those of typical or 
giant ellipticals obtained from optical stellar absorption 
hues (e.g. Ca H&K, Mg lb, CaT). This was based on 
the assumption that the near-IR stellar absorption lines 
or stellar band-heads probe the global properties of the 
ULIRGs. The results presented here and in Paper I 
imply this is not the case. The second is comparing 
the values of A/Dyn, computed from the CO CTo and 
near-IR photometry, with some fiducial stellar mass 
representative of the stellar mass function of galaxies. 
The earlier studies above all used m* — 7.07x10^" 
Mq (or 1.25x1 ^"'^ M (7) for the cosmology used here) 
from iCole et al.l ()2001h . In other words, because in the 
near-IR MDyn < n^*, ULIRGs cannot form gEs, let alone 
an average elliptical. They must be the progenitors of 
low-intermediate mass ellipticals. However, the value 
of m* used for com parison is a. c tually the larger of two 
possible values from iCole et all (|200lD . The other is m* 
= 3.43x10^° h-^ Mq (6.0x10^° M© for the cosmology 
used here). The larger value comes from using a Salpeter 



IMF which over-predicts the amount of low-mass stars, 
rather than a Kennicutt or "diet" Salpeter I MF which 
comp ensates for this e ffect ( see Section 6.1 of iBell et al.l 
[2001 iBeh fc de Jcmj [Ml . A variety of methods 
have converged towards m* ~ 3x10^° Mq (the value 
used here) which also appears to be the transition 
region between the blue cloud and red sequence and 
the threshold ab ove which AGN activ i ty is m ore likely 
to be found (e.g. IBell fc de Jonel [2Q0l[ : liBell et al. 2001 
iKauffman n et al.l 120031: iBaldrv et'Iil 120041 : iBeh et all 
2007; Bald rv et al.l 120081 ) In other words, the claim 
that (U)LIRGs form sub-m* ellipticals is partly based 
on the selection of the larger of two possible m* values. 
In the bottom panel of Figure 5 the ULIRGs straddle 
the m* line. This remains the same when near-IR 
data fo r other advan ced ULI RGs from [ Genzel et al.l 
(pOOll ): iTacconi et all {2002); Dasvra et al.l ()2006[ ) are 
considered. Instead, the real question raised by these 
earlier results is why the observed CTo.co of ULIRGs 
and near-IR half-light radii were inconsistent with 
observations at other wavelengths, including molecular 
gas masses, star-formation rates, high- velocity outfiows, 
etc. Similar to the results for non-IR luminous mergers 
and LIRGs (e.g. RJ06a, Paper I), the kinematic and 
photometric properties of the ULIRGs measured at 
/-band are now statistically consistent with m* and 
larger ellipticals, including gEs. 

QSOs and ULIRGs have a great number of simi- 
larities, including: bolometric lumi nosities and space 
densities out at least z ~ 0.4 (|Soifer et al.l 119861 : 
iCanalizo fc Stocktonll2001[): an overla p in th e FIR-radio 
corre lation (e.g. iSanders fc Mirabell Il996t lYun et al 
200l: masses of ~ 10^-^° Mq, (e.g. ISanders' et al 



1989b; 'Chini et al."1997: 'Evans et al."2001: ScoviUc et aL 
2003 ; Bertram ct al. 2007; Evans 2009; Evans et al. 200St 
Xia et al.l [20121 ): post-starburst stellar populations in 



or near the nucleus along with tidal tails and peculiar 
morphologies indicative of a relatively recent gas- 
rich mer ging event ( e.g. 'MacKentv fc StocktonI 
Heckman et al.l 119861 : Hutchings & Ncff^ 119881 . 
Guvon et a l. 2006; Canalizo ct al. 2007; Bcnncrt et al 
20081 [Ramos Almeida et al. MIL; Tadhunt er et al 



and an overlap in the distribution of Lir (12.24 
± 0.44, 12.25 ±0.47, and 12.17 ± 0.16 for RLQs, 
RQQs, and ULIRGs, respectively), in which the Null 
Hypothesis cannot be rejected. However, the sizes 
and masses of ULIRGs and QSO host galaxies have 
previously been reported as significantly different. The 
/-band dynamical results presented here alleviate this 
discrepancy. There is now a much stronger dynamical 
link between ULIRGs and QSO host galaxies. The 
strongest result comes from the KR, in part, due to 
the large sample size of QSOs. The 2D KS test for the 
distribution of objects in the Log /?eff-<Mi>off plane, as 
well as the two-sided KS tests for each parameter rule 
out statistical differences between ULIRGs and QSOs 
(whether grouped together or compared separately as 
RLQs and RQQs). While past comparisons have relied 
on the KR to reject the notion that ULIRGs evolve into 
QSO host galaxies, the same comparison here at /-band 
strongly supports the S88 paradigm. 

However, it is still important to compare the kinemat- 
ics of ULIRGs and QSOs. Unfortunately, the tradeoff 
for doing so is the significantly smaller sample size of 
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QSOs, including the loss of comparing ULIRGs with 
only RQQs. The results for the FP are consistent with 
those of the KR. QSOs (either taken together or just 
RLQs) and ULIRGs show no statistical difference in 
their distribution in FP parameter space. It is only 
when the parameter CTo is considered alone that things 
become less clear. There is a weak statistical difference 
between RLQs and ULIRGs for this parameter (see 
Table CI), primarily due to the large RLQ errors. In 
the case of MDyn the differences are stronger and the 
RLQ fTo errors do not affect the results significantly 
(see Table CI). How can these results be reconciled 
with those from the KR (and its individual parameters) 
which uses the larger photometric QSO sample? Is 
cTo really different for RLQs and ULIRGs, or does 
it appear to be different because the 6 RLQs in the 
kinematic sub-sample happen to be non-representative 
of the larger RLQ sample? To test this, the 2D and 
standard two-sided KS tests were re-run for the KR and 
its parameters between the ULIRGs and the kinematic 
QSO sub-sample only. Just one difference emerges from 
the results listed in Table CI. The Null Hypothesis 
can be rejected for the Log Res comparison between 
ULIRGs and the 6 RLQs. This explains the rejection 
of the Null Hypothesis for Afnyn between RLQs and 
ULIRGs because (Tq and Log i?eff are used to compute 
the mass. Although the kinematic results for the RLQs 
are uncertain, when the QSO hosts are considered as 
a single population, it suggests kinematic similarities 
exist between them and ULIRGs. Overall, these results 
clearly demonstrate the need for more measurements 
of (To in RLQs and RQQs to confirm the results from 
the KR and probe whether the kinematic differences 
between RLQs and ULIRGs are real (and if any exist 
between ULIRGs and RQQs). 

Although beyond the scope of this paper, these results 
also raise a conundrum in regards to whether RQQs and 
RLQs are dynamically different. The Null Hypothesis 
can be rejected when their luminosities are compared. 
This implies that their host masses are different (assum- 
ing some M/ L transformation from stellar population 
models and that the ages of the stellar populations in 
RLQs and RQQs are the same). However, their 2D 
distributions in the KR are not statistically different 
(nor are the two parameters when each is compared 
separately). Since the KR is a projection of the FP, 
which in turn is related to the correlation between M 
and M/L, it implies that RLQs and RQQs are not 
dynamically different. 

7. SUMMARY & FUTURE WORK 

The main results of this paper are summarized below. 

1) The o"-Discrcpancy. first reported in RJ06a and 
quantified in Paper I for LIRGs is shown to extend to 
ULIRG luminosities. The (Tq measured from the CaT 
stellar absorption lines are systematically larger than 
those obtained from the near-IR CO band-heads. With 
the addition of ULIRGs the correlation between (TFrac 
and Log Ljr remains unchanged from Paper I. We posit 
that for the single nuclei (U)LIRGs presented here and 
in Paper I, this relationship results from feedback pro- 
cesses that cause monotonic aging and dimming of the 
YCSD population and the clearing out of the dusty, star- 
forming medium. 



2) At /-band, ULIRGs are nearly an order of magni- 
tude more massive than previously measured in the near- 
IR, and are consistent with ellipticals ranging from m* 
to gEs. All of the ULIRGs presented here lie closer to 
the Fundamental Plane and Kormendy Relation than in 
near-IR studies. 

3) At /-band, the M/ L values of ULIRGs appear to 
indicate the presence of an old, evolved stellar popula- 
tion, similar to quiescent ellipticals. Yet in the near-IR, 
ULIRGs reflect much younger populations, matching the 
well established observations of significant quantities of 
molecular gas and high rates of SFR. 

4) At /-band ULIRGs are dynamically similar to QSO 
host galaxies, further supporting the S88 paradigm. 
ULIRGs are statistically consistent with the positions of 
both RLQs and RQQs on the Fundamental Plane and 
Kormendy Relation. This result uses the same meth- 
ods of comparison which in the past have been used to 
demonstrate that ULIRGs do not evolve into QSO host 
galaxies. However, when the ULIRGs are compared with 
the kinematic sub-sample of RLQs, there is a statistical 
difference (tlw^ Null Hypothesis can be rejected) for ito, 
MDyn, and Log /?off. The impact of this difference is 
weakened by two caveats; the large errors in ctq for the 
RLQs which affects the KS tests (sec Table CI;) and 
the contradiction which arises from the Niill Hypothesis 
not being rejected for ULIRGs and RLQs when the full 
sample of 28 RLQs are considered. These results demon- 
strate the need for more CTo measurements for RLQs (to 
resolve the contradiction) and RQQs in order to make a 
viable kinematic comparison with ULIRGs and confirm 
the results presented here. 

5) Finally, an homogenized /-band sample of RLQ and 
RQQ host galaxies are presented here which can be used 
for future dynamical studies. The QSO hosts can be used 
as cither a control sample for further comparisons with 
IR-luminous systems or as a representative sample itself 
for future studies of QSO host galaxies. 

These are the first results from a much broader survey 
to establish an accurate mass distribution for ULIRGs 
and to rc-cvaluatc how these systems fit into the broader 
picture of the formation and evolution of elliptical galax- 
ies and QSOs. It is clear that more kinematic (stellar (Tq) 
observations are needed for RLQ and RQQ host galaxies 
in order to confirm the results presented here. However, 
the results presented so far are consistent with both the 
Toomre Hypothesis and QSO evolution scheme presented 
in S88. Future work will focus on a multi-wavelength 
approach which will contimie to use optical observations 
to measure MDyn and near-IR observations to probe the 
central kpc in ULIRGs. A comparison of central black 
hole masses (M.) will also be made, in order to deter- 
mine the location of ULIRGs with respect to the M.-rr 
relation. These will require observations to infer M, in- 
dependent of a or host luminosity. Assessing detailed 
and accurate properties of local ULIRGs is key to for- 
mulating a better understanding of similar systems (e.g. 
dust obscured galaxies, sub-millimeter galaxies, etc..) at 
higher redshifts. 
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Table 1 

ULIRG Galaxy Sample 



Galaxy 


R.A 




Dec. 


z 


Log Lia 


E{B - V) 


Name 


(J2000) 


(J2000) 






(mag) 


IRAS F02021-2103 


02 


04 


27 


-20 49 41 


0.116 


12. 02''''' 


0.020 


IRAS 05189-2524 


05 


21 


01 


-25 21 45 


0.043 


12. ll'' 


0.026 


IRAS F10378-I-1108 


10 


40 


29 


+10 53 18 


0.136 


12. 26='''' 


0.029 


IRAS 11387+4116 


11 


41 


22 


+40 59 51 


0.148 


12.03''''=''* 


0.015 


IRAS 12540+5708 


12 


56 


14 


+56 52 25 


0.042 


12. 48=^ 


0.009 


IRAS 17208-0014 


17 


23 


21 


-00 17 01 


0.043 


12. 34=^ 


0.304 


IRAS 19542+1110 


19 


54 


35 


+11 19 02 


0.064 


12. 03''''' 


0.199 


IRAS 23365+3604 


23 


39 


01 


+36 21 09 


0.064 


12.13='''' 


0.096 



Note. — R.A. is in units of hours, minutes, and seconds, and Dec is in 
units of degrees, arcminutes, and arcseconds. Ljr is the 8-1000 /.tm total flux 
measured from the 12, 25, 60, and 100 ^m IRAS passbands. (a) All or some 
IRAS fluxes from Sanders et al. (2003); (b) 12/xm flux from WISE; (c) All or 
some IRAS fluxes from iMoshir fc et al.i tl99a '): (d) 22/im flux from WISE. 



Table 2 

ULIRG Observation Log 



Galaxy 


Imaging 


Imaging 


ESI 


ESI 


ESI 


Name 


Camera/ 


Integration Time 


Integration Time 


RA. 


Slitwidth 




Filter 


(sec) 


(sec) 


(dcg) 


(arcsec) 


IRAS F02021-2103^ 


EST WFPC2/F814W^ 


800 


3600 


55.0 


1.0 


IRAS 05189-2524*^ 


i/5T ACS/F814W' 


730 


900 


0.0 


1.0 


IRAS F10378+1108'' 


SDSS/Sloan z 


54.1 


1800 


0.0 


1.0 


IRAS 11387+4116'' 


SDSS/Sloan z 


54.1 


1800 


0.0 


1.0 


IRAS 12540+5708'^ 


i/5T ACS/F814W' 

EST ACS/F814W',NIC2/F160W° 


830 


900 


45.0 


0.75 


IRAS 17208-0014=' 


720,224 


3900 


40.0 


1.0 


IRAS 19542+1110=^ 


F5T ACS/F814W'' 


720 


3600 


0.0 


1.0 


IRAS 23365+3604=^ 


EST ACS/F814W'',NIC2/F160W'' 


750,2496 


5400 


300.0 


1.0 



Note. — (a) ESI P.I. Rothberg; (b) ESI P.I. Sanders; (c) _ff5T Program 6346, P.I. Borne; (d) HST Program 10592, P.I. 
Evans; (e) EST Program 7219, P.I. Scoville; (f) EST Program 11235, Surace. The WFPC2 observations were centered 
on the WF3 chip. 



Table 3 

ULIRG Rest-frame /-band Properties 



Galaxy 
Name 


Mi 
(mag) 


-Rcff 

(kpc) 


</iI>cff 

(mag arcsec"^) 


CaT (To 
(km s-i) 


CaT Vq 
(km s~^) 


CaT 
Template Star 


(xlO"M0) 


IRAS F02021-2103 


-23.61±0-05 


951±0.23 


19 84+0.05 


209a±8 


34679±* 


G8III HD 35369 


5 79±0.46 


IRAS 05189-2524 


-22 71±0 09 


2.59±o.io 


17.93±oi2 


265='±'^ 


12869*^^ 


KOIII HD 206067 


2.54+0.15 


IRAS F10378+1108 


-22'96±010 


3.09±oi4 


18.05±°-23 


280''±" 


41007±9 


GlIII a Sge (HD 185758) 


3.38±0-3i 


IRAS 11387+4116 


-23.01±0-0i 


3 22±o.07 


18.09±0" 


198b ±9 


44575*^^ 


MOIII A Dra (HD 100029) 


1.76±oi6 


IRAS 12540+5708 


_23.42±0.10 


5.88±0-27 


18.98±°13 


34Qc±9 


12584±io 


G5II HD 36079 


9.86±0-68 


IRAS 17208-0014 


-22.88±''-" 


9 37±o.50 


20.54±0-i2 


2gia±5 


12798±"' 


KOIII HD 206067 


8.90±0.58 


IRAS 19542+1110 


-23.52±009 


o;77±o.03 


1443±o.i4 


IggaiS 


18718±5 


G8III HD 35369 


0.31±002 


IRAS 23365+3604 


_23.i6±o.08 


416±o.i9 


18.50±oio 


22ia±6 


19310±5 


K1.5III a Boo (HD 124897) 


2.83±o.i9 



Note. — Fluxes are in VEGA magnitudes. The values have also been c orrected for Galactic Reddening using dust maps from 
ISchlegel erall H1998I 1 and scaling from Table 6 in Schlaflv & FinkbeiniB l|20111 ). assuming _Rv = 3.1 (FitzDatrick 1999). The Ax scaling 
factors used were: Ap8i4-w (ACS/WFC) = 1.52; v4pgi4w (WFPC2) = 1.54; and = 1.26. Circular isophotcs were used for the 
F8 14 W photometry. Elliptical isophotes were transformed to equivalent radii for SDSS z-band photometry on IRAS F10378+1108 and 
IRAS 11387+4116. (a) Measured in 1.53 kpc diameter aperture; (b) Corrected to 1.53 kpc diameter; (c) Measured 2.08 kpc NW from 
the nucleus and corrected to a 1.53 kpc diameter. 
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Table 4 

Near-IR ULIRG Properties 



Galaxy 


CO o-o 


Moyn via CO 


R,B {F160W) 




Name 


(km s~^) 


(xlO"M0) 


(kpc) 


(mag) 


IRAS F02021-2103 (A) 


143±21a 


1.09 


3.85^^ 


-24.89' 


IRAS 05189-2524 (B) 


;^3;^±16a 


0.11 


0.48=* 


-23. 96* 


IRAS 12540+5708 (C) 




0.20 


l.OS'* 


-24.22^ 


IRAS 17208-0014 (D) 


223±i5c 


0.94±o.i4 


I 3g±0.10e 


_24.43±0.18 


IRAS 23365-1-3604 (E) 


143 ±i5c 


0.57±0-" 


2;oo±0.22c 


-25.25±0-4l 



Note. — The letters A-E correspond to objects in Figure 5. Fluxes are in VEGA 
magnitudes. The values have also been corrected for Galactic Reddening using dust 
maps from Schlcgcl ct al. (1998) and scaling of j4Figow = 0.51, assuming fiv = 3.1 
||F itzpatrick , 1999 ^ . All values of (To,co have been corrected to a common aperture 
of 1.53 kpc. Some values in the table d o not incl ude errors because they were not 
available from the source materials, (a) IDasyra ct al. (2006), MOIII template star 
use d (HP 99817) ; (b) Tacconi ct al. (2002) K5Ib template star used (HD 200576); 
fc) IGenzel eFall (12001) MOIII template star used (HD 99817); (d) PSF-subtracted 
Sersic fit (n = 4) (Veilleux et al. 2006) from elliptical isophotes converted to equiv- 
alent radii using ellipticities derived from their Table 2; (e) Photomery measured 
using circular isophotes. ("f 1 iVeilleiix et al.l ((2006 ) . The values of (To, co ™ Column 
2 and Moy^ values computed in Column 3 are only for comparison purposes with 
the CaT derived properties in Table 3 and do not represent the global properties of 
the systems. 
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Table 5 

Kormendy Relation Fits 



Sample 


N 


Slope 


Intercept 


r.m.s. 


r 




Rest-frame /-band 


SDSS Ellipticals 


9255 




17.56±0-0i 


0.S5 


0.71 


SDSS Ellipticals 85 < CTo < 125 km s^^ 


1656 


2.76±o.oi 


17 75±o.oi 


0.S8 


0.76 


SDSS Ellipticals 125 < CTo < 165 km s"! 


S301 


2.69±o.oi 


17 SS*''-'^-'- 


0.S6 


0.82 


SDSS Ellipticals 165 < CTo < 225 km s'^ 


S524 


2.72±0-0i 


17^02±0-0i 


0.S5 


0.85 


SDSS Ellipticals 225 < CTo < 420 km s"! 


774 


2 37±0.02 


17.00±0-0i 


0.21 


0.87 


ULIRGs^ 


8 


4;50±o.i4 


15.68±0-09 


0.22 


0.98 


RLQs'' 


28 


3 5]^±0.52 


15.71±0-« 


0.54 


0.79 


RQQs^ 


25 


s'40±0-*4 


16.ii±o.33 


0.56 


0.84 


All QSO Host Galaxies'^ 


53 


S.25±0-32 


16.07±0-26 


0.57 


0.81 


Restframe /-band for Previously Published Samples'^ 


RLQs (O'Dowd et al. 2002) 


16 


3gl±0.66d 


15.S7±0-5S 


0.46 


0.83 


RLQs (D03) 


10 


3gg±0.47e 


15 77±o.42 


0.23 


0.93 


RLQs (DOS, excluding two outliers*) 


8 


3 05±o.30g 


16.42±0-29 


0.12 


0.97 


RQQs (DOS) 


9 


340±o.73h 


16.4S±0-58 


0.S7 


0.86 


RLQs (F04) 


7 


S.68±0-io' 


15.8S±0-08 


0.S7 


0.94 


RQQs (F04) 


7 


3 j^3±0.20i 


15.86±0-i^ 


0.45 


0.91 


RLQs + RQQs Combined (F04) 


14 


3;2o±o.07j 


16,oo±oo6 


0.48 


0.91 



Note. — As a comparis on, the earlier /?-band and l/-band KRs had slopes of 3.02 
HKormendvl 1197^ ) and 2.94 l IHamabe fc Kormendvl 119871 ). repsectively, using the surface 
brightness at the effective radius. Two fits each were made to the RLQ, RQQ, and 
RLQ-I-RQQ data to account for the absence of errors in the DOS sample. If errors were 
available a DWLSQ fit, which takes into account errors in X and Y, was used. Otherwise, 
the data points were equally weighted and a standard least-squares fit was used. The fits 
plotted in Figure 6 include all data for the ULIRGs, RLQs, and RQQs. (a) IRAS 19542-f 1110 
is considered an outlier because it is 6.3a- from the KR plotted in Figure 6, while all other 
ULIRGs are < 2a from the line. Excluding it changes the fit to: </j>cff 3.85^"'^^ xLog R^g 
+ 16.19±'' ", r.m.s.= 0.24 and r =0.96; (b) No errors are available for the DOS sample. The 
KR was fit using a simplified least-squares fit with equal weighting of errors for the depen- 
dent data points. When the DOS data are excluded and a DWLSQ method was used the fits 
change to: RLQs {N = 22): </i>eff 3.58±'' 0'^ xLog R^ff + 15.66±° '"^, r.m.s.= 0.57 and r 
=0.81; RQQs {N = 20): </i>cff S.79±0 ''9 xLog R^a + 15.15±° '"^, r.m.s.= 0.83 and r=0.83; 
All QSOs {N = 42): <ti>cfi 4.2S±°-''^ xLog R^g + 14.85±° °3, r.m.s.= 0.74 and r =0.82; 
(c) All of the QSOs used to derive the preivously published KRs are also part of the QSO 
Comparison Sample in Appendix B. The previously published fits used the surface brightness 
a4the effective radius, not the surface brightness within the effective radius (</i>ofr), as used 
here. The original fits used different cosmologies, fc-corrections, transformations to different 
rest-frame filters, and different fitting techniques. They have been refit here to provide a clean 
comparison; (d) Original published restframe /?-band published slope = 2.75 , (e) Origi- 
nal published observed F675W-hand pubUshed slope = 3.98±'^ ''i; (f) PC 1004+130 and OX 
169; (g) Original published restframe /?-band published slope = 3.19^'''®^; (h) F675W-hand 
published slope = 2.99^"'^*; (i) RLQ and RQQ samples from F04 were not fit separately; 
(j) Original published restframe V-band published slope = S.SS^'^'^. 
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Figure 1. Comparison of the central velocity dispersions (o"o) 
measured from the Calcium II Triplet (CaT) stellar absorption 
lines at A ~ 0.85/im (x-axis) and those measured from the CO 
band-head at 1.6 or 2.3/im (y-axis) for 5 ULIRGs. 
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Figure 2. Values of trpraci which quantifies the relative difference 
between (Jo,CaT and (Jo, CO for each advanced merger, are compared 
with the independently measured quantity Log Lir. The horizon- 
tal dotted line represents (Jprac = and the diagonal dashed line 
shows the weighted least-squares fit to the data. Also shown is 
the Pearson Correlation Coefficient (r). Values of r range from 
-1 (strong anti-correlation) to +1 (strong correlation). There is a 
strong correlation between o-jy^c and Log Ljr. 



bJO 
O 




5.0 5.5 6.0 6.5 7.0 

Log Og (km s'^) + 0.22 <)-i.,>eff (mag arcsec"^) 

Figure 3. The diagonal solid line plot ted in this figure rep- 
resents the /-band Fundamental Plane (H yde Hi Berna rdi 200913) 
seen edge-on. It is transformed from the Sloan i-band filter to our 
preferred cosmology. Plotted in this figure are: advanced ULIRGs 
(filled circles), the host galaxies of Radio Loud (open squares), and 
Radio Quiet (open fiattened diamonds) QSOs, and the SDSS com- 
parison sample of 9,255 ellipticals (open circles in four colors). In 
this, and subsequent plots, the SDSS ellipticals are shown binned 
into four groups based on ao- 1) 85 km s~^ < o"o < 125 km s~^ 
(1656 blue); 2) 120 km s"! < (To < 165 km s"! (3301 green); 3) 
165 km s-l < (To < 225 km s'^ (3524 orange); and 4) 225 km 
s-l < (To < 420 km s-i (774 red). The ULIRGs and QSO host 
galaxies have values of cto which would place them in either the 
3rd or 4th bin. Error bars are plotted for the ULIRG sample and 
QSOs (where available), along with the median errors for the SDSS 
sample. Over-plotted are the 1, 2, and 3cr scatter of the FP (dotted 
diagonal lines in dark to light grey). 
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Figure 4. This figure compares (MDyn) and M/L values at I- 
band. Plotted here are the ULIRGs, SDSS comparison ellipticals 
and QSO host galaxies (same symbols used in Figure 3). The over- 
plotted solid vectors in the panel show the temporal evolution of 
M/L for a burst single stellar population (SSP) with solar metal- 
licity and a Kroupa (solid black line) or Salpeter (grey line) IMF 
(M05). The value of M for the SSPs represents the stellar mass 
(m*) while the masses plotted for the galaxies are the total virial 
dynamical mass (MDyn)- The horizontal placement of the M/L 
vectors are for display only. The vertical dotted line represents the 
mass of an m* galaxy (^ SxlQi" Mq). Error bars are plotted for 
the ULIRG sample and QSOs (where available), along with the 
median errors for the SDSS sample. 
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Figure 5. Shown here is a plot similar to Figure 4, comparing 
the same set of ULIRGs at /-band (top) with their parameters 
at F160W {bottom). The top panel uses /-band photometry along 
with o"o,CaT and the bottom panel uses HST/Af/CS F/ 60 photom- 
etry with (To CO • The 5 ULIRGs plotted in both panels arc: A = 
IRAS F0202i-2103, B = IRAS 05189-2524, C = IRAS 12540-f 5708, 
D = IRAS 17208-0014, and E = IRAS 23365-1-3604. The solid vec- 
tors in both panels represent the temporal evolution of M/L at 
/-band {top) and F160W {bottom) for a SSP with solar metallicity 
and a Kroupa (solid black line) or Salpeter (grey line) IMF (M05). 
Once again, the Mfrom the SSPs represents the stellar mass (m*), 
while for the plotted galaxies it represents MDyn- The horizontal 
placement of the M/ L vectors in both panels are for display only. 
The vertical dotted line in both panels represents an m* galaxy 
(~ 3x10^° Mq). In the bottom panel ULIRGs A, B, and C have 
been PSF-subtracted (FWHM = 0.14 arcsec pixel"^ correspond- 
ing to 0.27, 0.11, and 0.10 kpc for A, B, and C, respectively). This 
has likely removed some of the star-formation contribution to the 
galaxy luminosity. Therefore, their vertical positions in the bot- 
tom panel are upper limits suggesting their ages and M/ L may be 
smaller than what is shown. 



APPENDIX 

A. SLOAN DIGITAL SKY SURVEY 
COMPARISON SAMPLE SELECTION 
CRITERIA 

The selection criteria used to generate the SDSS 
comparison sample are descri bed below. They ar e 
based upon the criteria used in iBernardi et al] (|2003b[ ): 
iHvde &: Bcrnardi. (,2009a, b) to select early-type galaxies 




0.5 1.0 
Log Rrf, (kpc) 



0.5 1.0 

Log Rrff (kpc) 



Figure 6. Shown in both panels is the Kormendy Relation (KR), 
a photometric projection of the Fundamental Plane. The left 
panel compares photometric properties of the sample of 8 advanced 
ULIRGs with a larger sample of 28 RLQ and 25 RQQ host galaxies. 
The diagonal solid black line plotted in this figure represents the 
/-band KR computed from the SDSS comparison sample of 9,255 
Ellipticals {r.m.s. = 0.39 dex). Error bars are plotted for the 
ULIRG sample and QSOs (where available), along with the me- 
dian errors for the SDSS sample. The diagonal dashed line shows 
the range of </xi> and Log R^g for an L* galaxy of constant lu- 
minosity. The dotted grey lines (from dark to light) represent 1, 
2, and 3xr.m.s. of the KR. The r.m.s. of the various samples 
are: ULIRGs = 1.01 (0.52 without IRAS 19542-f 1110), RLQs = 
0.95, RQQs = 0.91, all QSOs = 0.93. Plotted in the right panel 
are the derived fits to the ULIRGs, RLQs, and RQQs from Table 
5. The filled regions for each sample display the fit including the 
itlcr errors in the slope and intercept. The diagonal solid black 
and dotted grey lines represent the KR and the 1-3 X r.m.s. from 
the fit. 



(including SO or lenticulars) for deriving the Fundamen- 
tal Plane in several bandpasses. The criteria are divided 
into photometric (items 1-7) and spectroscopic (items 8- 
12) categories. The selection criteria required that the 
elliptical galaxy must be present in both the photometric 
(PhotoDbjAll) and spectroscopic (SpecDbj) databases. 
Only Sloan i-band photometry was used for the photo- 
metric catalog. No magnitude or flux limitations were 
imposed on the sample. 

The biggest differences bet ween the select io n crite - 
ria used h ere and those in IBernardi et al.l (|2003bD : 
IHvde fc Be rnardi (2009c3|B) are: redshift (z < 0.15 here 
vs. 0.36); restricting the sample to elliptical morpholo- 
gies only (affected by InlDev and InlExp; and ti > 85 km 
s~^ vs 60 km s~^ to avoid t he effects of instrument al res- 
olution (see Appendix B in IBernardi et al.l (|2003bf )'l. No 
parameters or values for those parameters were selected 
that would exclude any elliptical in the comparison sam- 
ple from being a part of the IIB09 sample. 

1) MODE — 1 selects the primary photometry for each 
object, rejecting possible duplications as wells as objects 
flagged as saturated, near the edge of a CCD. 

2) PARENTID = and nCHILD = rejects galaxies 
blended with other objects, or the child of a de-blended 
set of objects. Given the spatial resolution of the pho- 
tometry and spectral resolution of the data, de-blended 
superimposed objects may yield unrealistic measured pa- 
rameters. Not setting these parameters resulted in a 
number of double counted galaxies in which the same 
a is assigned to both parent and child galaxy, producing 
spurious results. 

3) FRACDEV = 1: defined as the fracti on of total flux 
which contributes to a de Vaucouleurs ()de VaucouleursI 
Il953f ) or r^/^ fit to the galaxy light profile. 

4) devAB > 0.6: defined as the axis ratio of the mi- 
nor to major axis for a de Vaucouleurs fit. Analysis 
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by iHvde fc Bernardil (|2009ai rbl) demonstrated that while 
FRACDEV = 1 ehminates most disk galaxies from SDSS 
photometric catalogs, a non-trivial amount remain. Val- 
ues of devAB > 0.6 significantly improve the removal of 
late- type galaxies. 

5) Type = 3: Morphological classification as a galaxy. 

6) InlExp & InlDev > -9999: InlExp and InlDev 
are maximum likelihood functions which estimates the 
best-fit model parameters convolved with an estimate of 
the seeing. Smaller values indicate a larger likelihood. 
The criteria cutoff assure that the measured values have 
meaning, as values of -9999 mean no data is available. 

7) InlDev at least 10% > InlExp: This criteria selects 
objects in which the likelihood of a de Vaucouleurs fit 
is 10% greater than an exponential fit. This is a recom- 
mended setting from SDSS for selecting elliptical galax- 
ies. 

8) specObj ID 7^ 0: reject objects without spectroscopic 
observations. 

9) SpecClass = 2: Spectral classification as a galaxy. 

10) eClass < 0: a one-dimensional classification 
of spe ctral type from Principal Component Analysis 
(|Yip et al. 2004) in which negative values correspond to 
absorption line galaxies with old stellar populations and 
positive values correspond to star-forming galaxies. 

11) snO, snl, and sn2 > 10.0: The values of a in DR7 
are measured by fitting the rest-frame wavelength range 
0.4-0.7 /im. The snO, snl, and sn2 criteria were selected 
to ensure a was measured from spectra with sufficient 
S/N. 

12) CT > 85 km s-i and < 420 km s^^: The SDSS 
spectra are re-sampled to a dispersion of log A = 10"^ 
dex pixel" ^ which corresponds to 69 km s^^. The actual 
spectral resolution varies from 85-105 km for galaxies 
in the SDSS spectra (Bernardi et al. 2003lJ) due to vari- 
ations as a function of wavelength. The DR7 (and DR6) 
a values differ from earlier data releases in that they 
aremeasured using a direct-fitting method (|Rix fc White! 
119921) with the assumption of a Gaussian profile, rather 
than a Fourier fitting routine. The latte r appears to bia s 
tr's < 150 km s~^ systematically higher (jBernar dil [2007h . 
The direct-fitting method is the same method used for 
measuring the Uo from the spectra of the ULIRGs pre- 
sented here although the profile shape is fit with a Gauss- 
Hermite polynomial rather than a Gaussian because the 
S/N is higher (see Section 3, RJ06a, and Paper I). The 
DR7 release notes also warn that a > 420 km s"^ are not 
reliable. To err on the side of caution we have selected 
the lower cutoff to be the approximate limit of, rather 
than below the instrumental resolution. 

Finally, for completeness, we apply a rest-frame radius 
correction for each elliptical galaxy in the SDSS compar- 
ison sample. Since early- type galaxies have color gradi- 
ents, yielding slightly larger Rp^ at bluer wave lengths, 
we use equation 6 from lHvde fc Bernardil ()2009a|) . which 
interpolates the observed radii from adjacent bands. The 
median correction is -1-0.012 kpc. 



B. DATA FOR THE COMPARISON SAMPLES 
OF RLQ & RQQ HOST GALAXIES 
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Table 1 

QSO Comparison Sample 



Galaxy R.A. Dec. z Log Ljr E(B -VY Camera/ Mi _Rcft <Ml>cff 

Name (J2000) (J2000) (L©) (mag) Filter (mag) (kpc) (mag arcsec'^) 
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0, 


.043 


WF3/F791 10 


-24, 


.81±0' 


02 


12.52±0-13 


19, 


.23±0' 


05 


HB89 2215-037 


22 


17 


47 


-03 32 38 


0, 


.242 


12, 


37b,3,4 


0, 


.106 


PC /F702W* 


-23, 


.88±0' 


02 


g79±0.06 


18, 


.84±0' 


08 



Note. — PC = Planetary Camera; WF2 = Wide Field 2; WF3 = Wide Field 3; WF4 = Wide Field 4. All rest-frame /-band photometry 
is in Vega magnitudes and has been corrected for Galactic Reddening. The Ax scaling factors used were: ^f606W = 2.41; ^f675W = 2.52; 
ApY()2W = 1 .94; j4p79iw = 1.74; ^Fgi4w = 1-54. The ijgff listed in this table are equivalent radii, (a) Galactic reddening values from 
HSchlegel et a l. 1998); (b) Fluxes from NASA/IPAC Sca n Processing and Integr ation tool; (c) 12 Aim flux f r om WI SE; ( d) 22fim flux from WISE; 
(e) computed only from 12 and 2 2Mm WISE fluxes; (i) llmpev fc Neugebaueil ||T983); (g1 iSanders et al.l l|1989al ): (hi IHaas et"an l(200Cf l based 
on ISO photometry; (i) IHaas et al. (2003) based on ISO photometry; (j) IRAS Faint Source Catalog; (1) upper limits {Sxr.m.s.) used for 
/12; (2) upper limits (3xr.m.s.) used for 725; (3) upper limits (Sxr.m.s.) used for /eo; (4) upper limits (3xr.m.s.) used for /loo- The IRAS 
flux es for 3C59 were computed from 3cr upper limits only and noted in parenthese s above, (f) HST/ WFPC2 Photometry from (Flovc Tet al.l 
[2004, Floyd, private communication}; _(*)_HST/WFPC2 Photometry from ( Hamilto n et al.l l2002l. [ 20081 , Hamilton, private communication); (|) 
HST/WFPC2 Photometry from Ounlp p et al. (2003). Although i?-band values are listed in the paper, they are actually values for the F675W 
fllter and were not transformed to Cousins ij-band; (*) HST/WFPC2 Photometry from iBahcall et all ([19971 1. 



22 



ROTHBERG ET AL. 



Table 2 

QSO Host Kinematic Properties 



Galaxy 
Name 



Optical (To 
(km s"^) 



(xlO"M0) 



Radio Loud QSOs 



PKS 0736+01 


348 ±«3 


28±7.62 


PG 1226+023 


339 ±5« 


13.54±5i5 


PG 1302-102 


386 


17.13±713 


PKS 2135-14 


310 ±106 


13 78±io.50 


4C +31.63'' 


318 ±« 


9 92±3-39 


PG 2349-014=^ 


279 ±'54 


19'.78±1009 


Radio Quiet QSOs 


HB89 0054+144 


167 ±" 


1.76±0-26 


PG 1444+407 


313 ±22 


6.81±i-07 



Note. — All (To values are measured using 
the Ca II H&K stellar absorption line over 
the 0.385-0.42 ^m wavelength range. All val- 
ues of (To have been corrected to a common 
aperture of 1.53 kpc. The average radii of 
the extracted apertures are listed in Table 3 
ofpVolf & Sheinis (200^). (a) The average ra- 
dius is _usod for multiple extracted positions 
ifWorf' £~Sherniir2008» . 



C. IMAGES, LIGHT PROFILES, & SPECTRA 
OF THE ULIRG SAMPLE 



IRAS F02021-2103 



IRAS 17208-0014 




IRAS 05189-2524 




IRAS 19542+1110 



- ^^^^^^^ 



IRAS 12540+5708 



IRAF 23365+3604 




0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 0.5 1.0 15 2.0 2.5 
Radius^'^ (kpc^'^) Ratiius^'^ (kpc^'^) Radius^^* {kpc^'^) 

Figure 2. Plotted here are the surface brightness profiles in the 
F8U ly-band (filled circles) for the 6 ULIRGs observed with HST 
using ACS/WFC or WFPC2. Also shown are profiles for two 
ULIRGs observed with HST using the NIC2 camera at F160W 
(open circles). Icr standard errors are over-plotted on each point. 
The surface brightness profiles are measured out to a S/N=3 over 
the background. All profiles are measured using circular apertures. 
The plotted points are equally spaced and linear, corresponding to 
3 pixels for ACS, 2 pixels for WFPC2, and 2 pixels for NIC2 (a 
radius of 9 pixels was used for IRAS 12540+5708 to compensate 
for the saturated/masked central region). The light dashed line 
in each plot represents the best-fit de Vaucouleurs r^/* fit to all 
of the data. The Xv^ of the best fit is shown in each panel. The 
two profiles plotted in the panels for IRAS 17208-0014 and IRAS 
23365+3604 have not been shifted. They represent the actual val- 
ues. These two ULIRGs also show significantly redder (/ — H) 
colors at R < 1 kpc th an the value of 1.77-1.79 expe cted from a 
typical elliptical galaxy l|Frogel et al.llT979 : IPahrelll999l ^. consistent 
with the results from Figure 13 of Paper I. 




' 1^19542+1110 . 



IRAS F02021-2103 



IRAS 05189-2524 




0.855 0.86 0.865 

Rest Wavelength (|im) 



0.85 0.855 0.86 0.865 

Rest Wavelength (jim) 



IRASF1037K4-110fi 



IRAS 11387-f4116 




IRAS 125404-570S 



0.855 0.86 0.865 

Rest Wavelength {[im] 



lliAS 17208-0014 



0.855 0.86 0.865 

Rest Wavelength (|jm) 



Figure 1. HST F8I4W or SDSS z images for the 8 advanced 
ULIRGs. The images are presented in reverse grayscale with a 
logarithmic stretch. Overplotted on each image is a horizontal 
solid bar representing 5 kpc. In each image North is up. 



— 0.8 
i 0.7 



IRAS 19542+1110 




0.855 0.86 0.865 

Rest Wavelength ([im) 



0.855 0.86 0.865 

Rest Wavelength (iim) 



IRAS 23365+3604 




0.855 0.86 0.865 

Rest Wavelength ([im) 



0.855 0.86 0.865 

Rest Wavelength (jim) 



Figure 3. Shown here are the CaT spectra of the ULIRGs ob- 
served with ESI on Keck-2. The solid black lines show the ac- 
tual spectrum, the light grey lines show masked bad pixels and/or 
emission lines. The thick dashed line shows the best-fit convolved 
template. Also shown are the positions of known emission and 
absorption lines within the wavelength range (Note: this does not 
mean that all of these lines are detected.) 
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Table 1 

Summary of KS Tests: Can the Null Hypothesis be Rejected at the 95% level? 



Comparison 
Sample 


Fundamental 
Plane (2D) 


era 


MY,y^-M/L 

(2D) 


MDyn 


Kormendy^ 
Relation (2D) 


Log Res'" 


<Ml>eff'' 


Log LiR 


RLQs 


No (89%) 


Yes (26%) 


Yes^ (72%) 


Yes" (87%) 


No 


No'' 


No 


No 


RQQs 










No 


No 


No 


No 


All QSOs 


No (100%) 


No (97%) 


Yes (27%) 


No (100%) 


No 


No 


No 


No 


SDSS Ellipticals 


Yes'' 


Yes^ 


Yes^ 


Yes'' 


No8 


No' 


No 





Note. — This table presents a summary of the ID and 2D KS tests performed between the distributions of the ULIRGs and 
the various samples listed in Column 1. The term "2D" refers to the two-dimensional KS test. Otherwise the KS test results 
presented here are for the standard two-sided KS test which compares the distributions of two empirical samples. Due to the 
large errors associated with the RLQ Uo values, the KS tests also have been run using Monte Carlo simulations in which each 
RLQ (To is randomly assigned a value of (To ± Acto . Test results above which include the RLQ cTo parameter show a parentheses 
indicating what percent of 10,000 Monte Carlo simulations the result occurs. Although IRAS 19542-1-1110 is 3.9(r outlier from 
the FP, excluding it from the KS tests in columns 2-5 does not change the results presented above, (a) Comparison with larger 
photometric sample of QSOs (28 RLQs and 25 RQQs); (b) No for 165-225 km s"' and 225-420 km s"' bins; (c) The Null 
Hypothesis can also be rejecte(i at the 99% level; (d) The Null Hypothesis can also be rejected at the 99% level but only 13% 
of the time; (e) The Null Hypothesis can also be rejected at the 99% level but only 12% of the time; (f) No for the 225-420 
km s~' bins; (g) The Null Hypothesis can be rejected when IRAS 19542-1-1110 is excluded. The Null Hypothesis cannot be 
rejected for the 165-225 km s"' and 225-420 km s"' bins whether or not IRAS 19542-1-1110 is excluded; (h) If only the 6 
RLQs in the kinematic sub-sample arc compared with the ULIRGs, then the Null Hypothesis can be rejected; (i) When IRAS 
19542-1-1110 is excluded the Null Hypothesis can be rejected. 
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D. COSMIC RAY REJECTION ALGORITHM 
AND BAD PIXEL MASK 

Most HST/ ACS programs employ cither CR-SPLIT, 
(two separate exposures at the same pointing), or 
muhiplc (e.g. > 2) dithered positions to remove cosmic 
rays (CRs) and artifacts. The CR-SPLIT mode is most 
appropriate for programs where the absence of data 
in the gap wih not impact the science (small targets 
or point sources). MULTIDRIZZLE compares the two 
images and flags pixels that have changed significantly 
between the two exposures. Dithering will "fill in" the 
chip gap and allow for the recovery of information in 
the chip gap. With at least 3 dither positions the same 
technique for removing CRs in the chip gap can be 
used. Program 10592 used a two position dithering 
scheme (ACS-WFC-DITHER-LINE) that shifts the image 
5 pixels in X and 60 pixels in Y. This fills the chip 
gap but with data from one exposure only (each gap is 
filled by information from the other exposure). Thus, 
the final images contained significant CR hits and hot 
pixels in the chip gap. Because the ULIRGs were not 
centered on either chip, but centered in the ACS/WFC 
FOV (aperture WFC) the chip gap runs through the 
outer regions of the galaxy, impacting the science data. 
Because more objects were affected by this than in 
Paper I, an improved and more automated algorithm 
was developed to create a bad pixel mask. First, a zero 
level background image was created by identifying the 
median background flux levels of the multidrizzled final 
image and replacing pixels at or below these values 
with a value of zero (using IMREPLACE). Second, the 
zero-level background image was passed through a 
median filter using a 15x15 pixel filter box, creating a 
new filtered image. Third, the zero-level background 
image was divided by the median filtered image. In 
this divided image, all pixels with flux values above 
the maximum pixel value in the nucleus were set to a 
value of 1 (bad), the remaining pixels were to set to a 
value of (good). This pixel mask proved successful 
for identifying saturated stars, diffraction spikes, and 
elongated CRs in the gap and areas covered by only one 
pointing. This pixel mask was then combined with one 
created from pixels flagged bad by MULTIDRIZZLE, and 
a pixel mask created from the positions of foreground 
stars and background galaxies. 
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